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NOMOGRAPHS FOR EVALUATION OF TENSILE 
DATA OF RUBBER! 


By SHIGEO KASE 


ABSTRACT 


Nomographs for the evaluation of tensile strength of rubber vulcanizates are 
described. The use of such nomographs is illustrated. The calculations necessary 
for preparing the nomographs for use with five samples are tabulated. 


Use of a nomograph to evaluate tensile data quickly is desirable when 
results of measurements of tensile strength of rubber are distributed doubly 
exponentially (1) since this is complicated for routine work. Nomographs for 
this purpose are described here with exclusion of detailed discussions. They 
will be valuable for purchasing, production, quality-control, and research staffs. 

The minimum-variance and unbiased estimation of averages from at least 
five or six observations is recommended in presenting data of tensile strength 
of rubber (2). This estimation is based on the method of least squares and is 
performed by calculation. If there is no computer available, as is often the 
case with routine work, the calculation is troublesome. The nomographs which 
will be explained in what follows are constructed in order to reduce this work 
without the uncertainty or uneasiness which there may be in graphical 
analysis (3). 

CONSTRUCTION OF THE NOMOGRAPH 


If the tensile strengths are ranked from S; to S, in order of decreasing 
magnitude, the minimum-variance and unbiased estimators for averages are 
given by (average equals sum of products of individual values and the re- 
spective weights) 


i=1 


For m equal to 5 the values of w; are in Table I (2). The nomograph for 
estimation of the mode, the most pertinent average for evaluating tensile 
data of rubber (2), is constructed from the relation: 


S = 0.4189S, + 0.2463S, + 0.16765; + 0.1088S, + 0.0584S;. 


1Manuscript received July 18, 1955. 
Contribution from Togawa Rubber Manufactory, Sumiyoshi, Osaka, Japan. 
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TABLE I 


WEIGHTS WHEN 1 Is 5 IN Ea. [1] 








Average Wi We W3 





Mode 0.4189 0.2463 0.1676 0.1088 
Median 0.2345 0.2487 0.2154 0.1754 
Mean 0.1285 0.2500 0.2429 0.2137 





All nomographs are constructed in the same way and have n plotting lines 
for tensile strength S; (¢ = 1,2,...,m),m — 2 auxiliary linest; (7 = 1,2,..., 
n — 2), and one estimation line § for the average tensile strength. In Figs. 1-3 
which may be used to estimate the averages from five measurements, the 
plotting lines are light solid, the auxiliary lines are dotted, and the estimation 
line is heavy solid. Since the nomographs in Figs. 1-3 are specially constructed, 
each line has the same ordinate. The graduation on each line is arbitrary but 
linear. In general the auxiliary lines need not be graduated. 

This gives advantages to the nomographs. When the range of observations 
(or the difference between the greatest and the smallest values among the 
tensile strength) is small, accurate estimation can be performed by graduating 
the ordinate widely so that the range may cover the scale. When the range is 
large the estimation becomes possible by compressing the scale. If the abscissas 
corresponding to all lines whose distances are equal or proportional to those 
in Figs. 1-3 are transferred to vertically long section paper, the nomograph 
will be capable of repeated application by moving the ordinate slightly up- 


wards or downwards. For convenience of reproduction of these nomographs, 


the abscissas of all lines are in Table II. 
From Table II, the abscissas of three auxiliary lines (t;—/;) and estimation 
line (S) are equal on the nomographs for mode, median, and mean. Nomo- 


STRENGTH KGM./cM.? 


TENSILE 


i%&Se2te ts $ 


Fic. 1. Nomograph for estimating the mode from five measurements of tensile strength 
of rubber. 





: TENSILE STRENGTH 


STRENGTH KGM./cM.? 


TENSILE 


Fic. 2. Nomograph for estimating the median from five measurements of tensile strength 
of rubber. 

Fic. 3. Nomograph for estimating the mean from five measurements of tensile strength 
of rubber. 


graphs can be overlapped on one sheet of section paper, if one allows for the 
complication and possible confusion caused by the scatter of the lines except 
the S; line. 


TABLE II 


EXACT ABSCISSAS OF ALL LINES IN THE NOMOGRAPHS FOR 
ESTIMATING AVERAGES FROM FIVE MEASUREMENTS 








Average Si So S3 Ss Ss 


™ 





Mode 0 0.68 5.22 9.90 19.39 
Median 0 0.49 3.49 6.23 10.19 
Mean 0 0.38 2.81 5.16 8.32 


— i ee 
bo bo bd 
Croco 





PROCEDURE FOR USING THE NOMOGRAPHS 


A summary of the method of using the nomographs is as follows: 

1. Arrange the five measurements of tensile strength in order of decreasing 
size, and then rank them from S; to Ss. 

2. Plot these values as ordinates on the corresponding plotting lines located 
on a linear scale at abscissas shown in Table II. 

3. Draw a line connecting points S; and S», and settle the point of inter- 
section between this line and the first auxiliary line 4. 

4. Determine the intersection between the line connecting the point chosen 
on the auxiliary line ¢; and point S; and the second auxiliary line fs. 

5. Obtain the intersection on the third auxiliary line ts by connecting the 
point on the second auxiliary line tf, and S,. 

6. Read the ordinate of the intersection on the §-scale, which is fixed by 
connecting the point chosen on the third auxiliary line t; and Ss. The final 
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ordinates are the averages by the minimum-variance and unbiased estimators 
relationship, Eq. [1]. 

The use of the nomographs is illustrated by an example. The observations 
are the same as in the previous paper (3), that is, the tensile strengths are: 


Si = 194, Sp = 193, S3; = 192, Sy = 174, and S; = 166 kgm./cm?. 


The estimation of averages of these observations is shown in Figs. 1-3. Esti- 
mated values of averages agree extremely well with the values calculated from 
Eq. [1] and the coefficients in Table I. 


Appreciation is expressed to Dr. H. L. Williams, Polymer Corporation 
Limited, for assistance with the preparation of the manuscript. 
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A FREQUENCY STABLE TRANSISTOR AUDIO OSCILLATOR 
OF VERY SIMPLE DESIGN! 


By W. D. Epwarps 


ABSTRACT 


Attention is called to the excellent frequency stability of a very simple tran- 
sistor oscillator. The frequency stability with respect to supply voltage, tempera- 
ture, load, and transistor interchange has been investigated. 


Transistor audio oscillators are simple to construct but their frequency 
stability often leaves much to be desired. The inherent stability of a number 
of simple oscillators has been determined. Point contact and junction tran- 
sistors have been used. The most stable circuit tested in the range 200—1300 
cycles per second is described in the following paragraphs. 

Presently available transistors are low impedance devices and their use at 
low audio frequencies with parallel resonant circuits appreciably reduces the 
Q of the circuit unless the L/C ratio is kept small. However, the large physical 
size of suitably valued mica condensers used for stability is not always desira- 
ble. The series resonant circuit with its low impedance at resonance is more 
suited to this type of oscillator. 

The most stable series resonant oscillator of simple design which has been 
built uses a point contact transistor and is shown in Fig. 1. Without the 
inductor, L, it is seen that the circuit of Fig. 1(a) is that of a relaxation oscil- 
lator. Even with ZL in the circuit, relaxation oscillations take place at low 
frequencies for which 24./LC ~ CR,. The diode shown in Fig. 1(6) together 
with the emitter forward characteristics maintains the emitter potential close 
to base potential. Emitter voltages far beyond the cutoff point which are 
associated with relaxation oscillators are eliminated. The sinusoidal mode of 
operation is extended to lower frequencies and a great increase in the stability 
of the circuit results. The transistor maintains oscillations by supplying 
current pulses to the series LC circuit during the conduction half cycle. These 
current pulses are represented by the major voltage peaks shown in Fig. 2. 
The minor peaks show the voltage which exists across the series tuned circuit 
during the period in which the emitter is effectively non-conducting. The 
magnitude ratio of the two peaks depends on the Q of the series resonant 
circuit and to a lesser extent on the collector load resistor R.. The maximum 
frequency stability obtainable from a simple circuit was the prime interest 
and no attempt was made to regulate the oscillation amplitude. 

The variation of frequency with collector supply voltage has been investi- 
gated using the circuits shown in Figs. 1(@) and 1(b). The results are shown in 
Figs. 3, 4, and 5. Without diode stabilization, Fig. 1(a), the circuit is most 

1Manuscript received July 21, 1956. 


Contribution from the Division of Radio and Electrical Engineering, National Research Council, 
Ottawa, Canada. Issued as N.R.C. No. 3726. 
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FIG, 2 


Fic. 1. Series resonant audio oscillator circuits: a. Conventional emitter current bias. 
b. Diode stabilization of oscillator frequency. 

Fic. 2. The voltage wave form observed at the collector with respect to ground. Circuit 
of Fig. 1b. L =7 henries. 


stable when large emitter bias currents are used. However, as expected and 
as seen from Figs. 3, 4, and 5, the circuit now operates only over a reduced 
range of collector supply voltage. The increase in frequency stability which 
occurs when the emitter bias source is replaced by a diode is clearly seen. 
With diode stabilization and a frequency of 1036 cycles per second a supply 
voltage change from 5 to 45 shifts the frequency by barely 0.5%. Similar 
results were obtained with a second Type 1768 transistor. From Table I it is 
seen that the oscillation frequency is reasonably independent of the transistor 


TABLE I 


A TABLE TO SHOW THE EFFECT OF TRANSISTOR INTERCHANGE ON THE OSCILLATOR 
FREQUENCY. CIRCUIT AS SHOWN IN Fic. 1) 


Transistor 


1768 
1768 
1768 
1768 
1768 
1768 
1768 
1768 


Frequency 1 in c.p.s. 


Frequency 2 in c.p.s. 


v 


1036.4 
1036.0 
1036.1 
1036.2 
1036.0 
1037.7 
1035.9 
1036.0 


i) 


crc 
NANO ONNN 
ST Ore ss 0 
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Nr orn wor 


SAS 
00. = | 


%@ 


toro 


bo DO DD DD ND ND ND DO 
. cHore 


Co 


1729 . 9, 1037.7 
1729 No. ; 59. 1036.7 
1729 Jo. 57. 1035.8 
1729 a 58. 1036.3 
1729 » % 57.S 1036.5 
1729 ; 58.¢ 1035.9 
1729 ce 58. 1036.0 
1729 5 9. 1036.4 
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COLLECTOR SUPPLY VOLTAGE 


Fics. 3a and 3b. Frequency versus collector supply voltage curves to illustrate the effective- 
ness of diode frequency stabilization. Nominal frequency: 250 c.p.s. Transistor T61. 
3a, R, = 33 K. 3b, R, = 3.3 K. 
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COLLECTOR SUPPLY VOLTAGE 


Fics. 4a and 4b. Frequency versus collector supply voltage curves to illustrate the 
effectiveness of diode frequency stabilization. Nominal frequency: 1000 c.p.s. Transistor T61. 
4a, R, = 33 K. 4b, R. = 3.3 K. 
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used. The maximum frequency shifts observed, and these with two types of 
transistor, are 0.2% and 1% at 1036 and 257 cycles per second respectively. 
Except where otherwise stated a collector supply voltage of 22.5, a load resist- 
ance of 6.8 K, and an inductance of 7 henries was used, both for these and for 
later results. 

In the initial determination of the frequency versus temperature relationship 
the temperature of the transistor alone was varied with the rest of the circuit 
remaining at room temperature. At a frequency of 235 cycles per second a 
change in frequency of approximately one per cent occurred as the temperature 
changed from 5° to 35°C. This variation was traced to the inductor L, a 
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COLLECTOR SUPPLY VOLTAGE 


Fics. 5a and 5b. Frequency versus collector supply voltage curves to illustrate the 
effectiveness of diode frequency stabilization. Nominal frequency: 10,000 c.p.s., L = 0.1 henry. 
Transistor T61. 

5a, R, = 33 K. 5b, R. = 3.3K. 
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10-henry filter choke having silicon-steel laminations. The transistor gain 
varied appreciably with temperature and this changed the current fed back 
through L from the collector to the emitter circuit. Because of the non-linear 
permeability, load current relationship of the inductor, this current variation 
changed the effective L and so the resonant frequency. This 10-henry inductor 
was replaced by a variable inductor (General Radio Type 1490B) which 
had much more linear characteristics. The 1490B is physically large and was 
not subjected to the temperature tests described below as was the rest of the 
circuit. The transistor was immersed in a silicone oil temperature bath and 
the results shown in Figs. 6(@) and 6(b) were obtained. In these figures the 
stabilization achieved by the diode is emphasized by comparing the results 


°Ro=2-7K | circuit 1(b) 
x Rco=6-8K 


Mm 





CIRCUIT I(a) 
Ie =0.I ma, 
Ro = 6.8K, 

Re =33K 
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TEMPERATURE © — 


Fic. 6a. Variation of oscillator frequency with temperature, with and without stabilizing 
diode. Transistor T55. Nominal frequency = 250 c.p 
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oRL= 68K) cigcuit i(b) 
XR, = 2.7K 


FREQUENCY IN CYCLES/SEC. 


CIRCUIT I(a) 
Ie =0.! ma, Ro = 68K, 


PH 


TEMPERATURE ‘C —> 


Fic. 6b. Variation of oscillator frequency with temperature, with and without stabilizing 
diode. Transistor T55. Nominal frequency = 1000 c.p.s. 


with those obtained for an emitter bias current of 0.1 milliampere. Between 
—60° C. and +45° C. the frequency is stable to within 0.1% at 1038 cycles 
per second, and to within 0.7% at 258 cycles per second. The latter figure may 
be reduced to 0.15% at 260 cycles per second by reducing the collector load 
from 6.8 K to 2.7 K. However, this increases the battery drain and reduces the 
starting reliability of the oscillator at the high frequencies, of the order of 
1300 cycles per second (L = 7h.). The temperature coefficient of the inductor 
will normally give rise to frequency variations in excess of these figures, and 
recognized methods may be used to compensate for these changes. The 
oscillation amplitude is approximately constant below 35° C. but decreases 
rapidly as the temperature rises above this point, and at 61°C. was about 
10 per cent of the amplitude at 20° C. The oscillator ceases to function within 
a few degrees of the highest temperature recorded. Very similar results were 
obtained with a second transistor, T97, Type 1768. 
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The long term stability of the oscillator appears to be excellent. At a fre- 
quency of 1038 cycles per second, and with a collector supply of 22 volts, the 
frequency did not change more than +1.5 parts in 20,000 over a 14-day period. 

The voltage wave form at the collector with respect to ground, obtained 
with the circuit of Fig. 1(0), is of the form shown in Fig. 2. This wave form is 
suitable for many purposes and the oscillator frequency is reasonably inde- 
pendent of loading at this point. A CR load in which C, and Ry, were varied 
over the range 0.01 uf. to 1 wf. and 5 K to 1 meg, respectively, was applied 
to the collector as indicated in Fig. 1(0). The maximum load condition of 
R= 5K and C, = 1 uf. caused maximum frequency shifts with respect to 
the unloaded circuit frequency of 0.6% and 0.15% for 250-cycle and 1000-cycle 
signals respectively. A large amplitude sine wave is available at point X, 
Fig. 1(), owing to the resonant rise across the condenser. However, this 
places the load in parallel with the reactive elements of the resonant circuit 
and will render the oscillator frequency dependent on the load. By restricting 
the loop gain, the oscillator (Type la) may be operated to give a low audio- 
frequency sine wave output at the collector, but the circuit adjustment is 
critical and this leads to instability. 

The measurements made on this oscillator are given in detail, for it is felt 
that the degree of stability obtainable with such a simple transistor circuit is 
not generally recognized. A very simple laboratory oscillator may be built 
around the circuit of Fig. 1(), and will be quite stable enough for most test 
purposes. By sacrificing frequency stability the collector signal amplitude 
for a given supply voltage may be considerably increased by replacing the 
load resistor R, with an inductor or with a transformer if desired. The circuit 
described is most suited to frequencies of less than about 1300 cycles per 
second. If greater frequency stability is desired the circuit will necessarily be 
more complicated, the first step being to stabilize the oscillator amplitude. 





ESTIMATION OF THE HEAT CAPACITIES OF ORGANIC LIQUIDS! 


By A. I. JOHNSON AND CHEN-JUNG HUANG 


ABSTRACT 


A method of estimating heat capacities of organic liquids at 20°C. based on 
the addition of group heat capacities is suggested. The group contributions were 
estimated for the most part from the data summarized by Chow and Bright. 
The method compares favorably with the correlation of Chow and Bright, 
being easier to apply and only slightly less accurate. The accuracy of the proposed 
method is much greater than that possible by assuming the atomic heat capacities 
suggested by Wenner. 


In recent years chemical engineers have become more aware of the need for 
rapid methods of estimating physical properties of substances for design 
calculations. Liquid heat capacity is one of these properties for which experi- 
mental data on new compounds is often lacking. 

Recently Chow and Bright summarized the existing methods of estimating 
liquid heat capacity (1); they suggested a new correlation involving the para- 
chor and molar refraction. Although their correlation works well (usually 
within 5%), the calculations are rather lengthy and the strain constants 
require some experience for proper use. A greater difficulty with the correlation 
is that the proposed equation involves two constants which apply toa particular 
type of compound. If the functional groups of an organic compound would 


place it in several groups at the same time, the heat capacities should be 
calculated using the several constants and presumably an average taken. 

Chow and Bright refer to attempts to calculate liquid heat capacity by 
statistical mechanical methods. Recently (4) Sakiadis and Coates reported 
success in using velocity of sound and infrared and Raman spectra data. These 
data may be difficult to obtain for many compounds, particularly for new 
compounds. 


A MODIFIED KOPP’S RULE 


Hougen and Watson (3) recommend the use of Kopp’s rule with Wenner’s 
atomic heat capacities (5). Chow and Bright point out that this ignores the 
constitutive properties of a molecule. The purpose of this paper is to report on 
an attempt to calculate group heat capacities which are additive and which 
reflect the constitution of the molecule. The authors have not been able to 
find any report of such a method, other than that of Wenner, in the literature. 

Molar heat capacities at 20°C. were calculated from the data summarized 
by Chow and Bright. 

Using homologous series of hydrocarbons, alcohols, acids, and esters, the 
contribution of the —CH.— group was found to be 6.32 cal./°C. (average of 
23 values). 

‘Manuscript received July 8, 19565. 
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From the homologous series of normal hydrocarbons (seven values) the 
contribution of the CH;— group was found to be 9.87 cal./°C. 

These values were then used to estimate the contributions of other groups 
shown in Table I. It should be emphasized that experimental values at 20°C. 
are scarce and that some of these values should be revised as more data become 
available. The values of —CN, —S—, and Br— were based on specific heat 
data taken from the Handbook of Chemistry and Physics (2). The list of 
benzene compounds of Chow and Bright was supplemented by five handbook 
values, resulting in a slightly higher group contribution for CsHs— than cal- 
culated from the data of Chow and Bright only. 


TABLE I 
GROUP HEAT CAPACITIES 


—-OH 
—-NHe (amines) 





| 


(esters) . | —Cl 
| —NO, 
o— —O— (ethers) 





fo 
en (ketones) 4.7 —Br 





—CN 3.9 | CH — 





USE OF THE SUGGESTED GROUP CONTRIBUTIONS 
The group contributions of Table I are used in the same manner as atomic 
contributions are used by Kopp’s rule, as illustrated by the following sample 
calculations: 


Isobutyl acetate 
O 
go CH, 
cH—ce OM 
Oot 
‘CH.CHs; 
3CH;— 3 x 9.87 29.61 
- 
tp 14.52 
—Chy 6.32 


Pa 
7‘. saa 
55.87 


The observed value is 53.27 cal./°C. per mole. 





Anisole 


C.H;— 
—O— (ether) 
CH;— 
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CsH;—O—CH; 
30.51 
8.37 
9.87 


48.75 


The observed value is 52.23 cal./°C. per mole. 


TABLE II 
COMPARISON OF CALCULATED AND OBSERVED HEAT CAPACITIES (20°C.) 





Compound 





Heat : 
capacity, 

gm-cal. 
(°C.) (gm.) 


M 


olal heat 


capacity, 


gm-cal. 


(°C.) (gm. mole) 
Obs. 


Calc. 








Normal hydrocarbons 
n-Butane 
n-Pentane 
n-Hexane 
n-Heptane 
n-Octane 
n-Nonane 
n-Decane 


Acids 
Formic acid 
Acetic acid 
Propionic acid 
n-Butyric acid 


Alcohols 
Methanol 
Ethanol 
n-Propanol 
n-Butanol 
Isopropanol 
Isobutanol 
Isoamyl] alcohol 


Esters 
Methyl formate 
Ethyl formate 
n-Propyl formate 
n-Butyl formate 
Methyl acetate 
Ethyl acetate 
n-Propyl acetate 
Methyl propionate 
Ethyl propionate 
n-Propy! propionate 
n-Butyl propionate 
Methyl butyrate 
Ethyl butyrate 
n-Propyl butyrate 
n-Butyl butyrate 


Ketones 
Dimethyl ketone 
Methyl ethyl ketone 
Methyl butyl ketone 
Diethyl ketone 
Dipropyl ketone 


0.549 
0.558 
0.527 
0.491 
0.505 

.503 

.493 


.512 
.488 
.473 
473 


.600 
569 
.571 

563 
.602 
.603 
.549 


07 
12 
‘ll 
ll 


.519 

508 
.459 
.459 
.468 
.457 
.459 
.459 
.457 
.459 
.459 
.459 
.457 
.459 
.459 


05 
4.08 
.10 
13 
.08 
.10 
13 
.10 
13 
5.16 
18 
13 
>. 16 
18 
21 


514 
.549 
553 
555 
0.550 


.08 
72.10 
100.16 
86.13 
114.18 


32. 
40. 
45. 
49. 
57. 
64. 
70. 


m Ors bo Go bo 


23.6 
29.; 
35. 
41. 


19. 
26. 
34.: 
41. 


36 


44. 
48. 


36 
37 


46 
40 
46 
53 


59. 
46. ¢ 
53. 
59. 
66. 


40. 
46.$ 
34. 
40.: 


Cor Ob 


Cow OW) 


Noe 


29.§ 
39.6 
55. 
47. 


ae. 


32. 
38. 
45. 
51. 
57. 
64. 
70. 


mre COON OC 


22. 
29. 
35. 
41. 


NWwWOs 


20. 
27. 
33. 
39. 
36. 
42. 
48. 


Cun ocn co 


27.§ 
34.: 
40. 
46.§ 
34.: 
40 
46. 
40. 
46.$ 
53. 
59. 
46.§ 
53. 
59.£ 
65.§ 


34. 


47. 


Molal heat 
capacity 
Deviation, calculated 
% Chow 
and 


Bright 


32. 
39. 
44, 
51. 
58. 
65 
ca 


00 Go “IER OOK SI 


24. 
29. 
35. 
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19 
26. 
34. 
42. 
37. 
43. 
48. 
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TABLE II (Concl’d) 
COMPARISON OF CALCULATED AND OBSERVED HEAT CAPACITIES (20°C.) (Concl’d) 


Heat Molal heat 
capacity, Mol. capacity, 
Compound gm-cal. wt, gm-cal, Deviation, capacity 
(°C.) (gm.) (°C.)(gm. mole) 7 calculated 

RED Chow and 
Obs. Calc. Bright 


Nitriles 
Acetonitrile 542 41. 22.2 23. >. § 22. 
Propionitrile .538 E 29.6 30. 6 29. 
Butyronitrile 548 ; 37.9 36. 3. 36.8 
Valeronitrile .521 3. 43.3 42. x 44. 





Amines 
Methylamine .765 31. 23. 25. 5. 24.; 
Ethylamine 691 31. 31.4 31. 
Propylamine .663 39. 37. 3.8 40. 


Chlorides 
Methy! chloride .371 d 18. 18.4 =i 18. 
Ethyl chloride . 369 ‘ 23. 24. ; 25.6 
Propyl chloride 404 78. dl. 31. F “ 
Butyl chloride 452 92.£ 41. 37. 6 40. 
Methylene dichloride . 287 84.$ 24. 23.£ 3. 24. 
Chloroform . 226 119.36 aT; 31. 5. 4 ‘ 


Esters of chloroacetic acids 
Methyl chloroacetate 0.383 108. 5¢ 41. 
Ethyl chloroacetate 0.399 122.£ 48. 
Propy! chloroacetate 0.415 136. 
Methyl dichloroacetate 0.311 143. 44. 
Ethyl dichloroacetate 0.330 157. 
Propyl dichloroacetate 0.342 171. 
Methyl trichloroacetate 0.268 177.4 47. 
Ethyl trichloroacetate 0.286 191. 
Propyl trichloroacetate 0.298 205. 4 61. 


39.¢ 
45.6 
61.{ 
47. 
53.¢ 
59.6 
52.§ 
59. 
65. 


or 


_ 
SIO RR NOITU 


bo 0 & Cn 00 


Benzene compounds and others 
Toluene 0.364 92.1: 33. 
Ethyl benzene 0.402 106.1! 42. 
Propyl benzene 0.400 120. 48 
Butyl benzene 0.428 134. 1! 57. 
Benzyl alcohol* (20-100°) 0.511 108. 1: 55. 
Benzonitrile (22—186°) 0.441 103. 45. 
Benzyl chloride (0°) .323 126. 40. 
Phenol (14-26°) .561 94. 52 
Benzophenone 383 182. 69 
Aniline 50 93. 46.! 
Allyl benzoate 388 162. 62.§ 
Allyl butyrate 451 128. 57. 
Allyl propionate 451 114. 51. 
Isoamyl propionate .459 144 66. 
Isobutyl acetate 459 116.06 53 
Diethyl malonate 433 160 69.: 
Anisole 483 108.1: 52. 
Diethyl ether 531 74 39 
Diethyl sulphide 477 90. 43 
Bromobenzene 231 157.02 36.5 
Nitromethane 412 61.04 25 
Nitrobenzene 339 123.11 41 45 9. 
Cyclohexanol 42 100.16 42. 48 14 
Benzene .406 78.11 31.7 32.5 2 


40. 
46. 
53. 
59. 
47. 
44. 
45. 
41. 
82. 
49. 
67. 
57. 
51. 
68. 
55. § 
61. 
48. 
40. 
43 
34. 
25 
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*The heat capacities of benzyl alcohol and following compounds are from the Handbook of 
Chemistry and Physics. 
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Table II compares calculated and observed molal heat capacities at 20°C. for 
84 organic liquid compounds. The average deviation was 5% with greatest 
deviations observed with benzene derivatives and esters of the chloroacetic 
acids. The values obtained by Chow and Bright are listed for comparison. 
Values calculated using the atomic contributions suggested by Wenner did 
not agree well with the observed heat capacities and are not shown. 

It is interesting to note fairly close agreement in the calculated value for 


benzene, based on six (CH groups. 
ie 


It should be noted that in formic acid and in the formates the contribution 
of hydrogen was taken as 3.55 (9.87 for CH;— minus 6.32 for CH2). 


CONCLUSIONS 


This paper proposes group heat capacities which may be used to estimate 
liquid heat capacities at 20°C. The method appears to be as useful as that of 
Chow and Bright in that it may be applied much more rapidly yet with suffi- 
cient accuracy for design purposes. 

Although comparison values are not shown in Table II, the suggested group 
contributions are more accurate than the atomic heat capacities of Wenner. 

No attempt was made to consider the temperature dependence of liquid 
heat capacity; Chow and Bright propose a method for extrapolating heat 
capacities to other temperatures which is probably satisfactory for design 
purposes. 
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THE CHEMICAL BLEACHING OF HYDROCARBON OILS! 


By J. R. FARNAND, A. E. McILHINNEY, W. S. PETERSON, 
P. E. GISHLER, AND I. E. PuDDINGTON 


ABSTRACT 


A method of bleaching hydrocarbon oils is described wherein the colored 
portions are selectively oxidized. In its most convenient form, air is blown 
through the oil at temperatures of 150-200°C. in the presence of metallic 
sodium. After a suitable period, the sludge formed is coagulated with water and is 
separated from the bleached oil by gravity settling or by centrifuging. The 
recovered oil has good oxidation stability and possesses detergent properties. 


INTRODUCTION 


The transparency of oleaginous materials is frequently improved by con- 
tacting them, at elevated temperatures, with adsorbents such as silica gel, 
alumina, fuller’s earth, bentonite, and bog iron. As an alternative procedure, a 
chemical method that is applicable to hydrocarbon oils has been examined 
recently and the results presented here indicate its possible utility. 

Under suitable conditions, it is possible to oxidize the colored components 
in oils selectively, and after removal of the oxidation products a material with 
a considerably improved color remains. Hydrogen peroxide, air, and oxygen 
at temperatures of 100-200°C. are suitable oxidizing agents while a number 
of alkali metals, their alloys or their salts, when added in small quantities 
improve the selectivity of the oxidation to a satisfactory degree. In some 
instances oils approaching medicinal oils in color, with only a slight bluish 
tint remaining, may be produced from either new unbleached raw material or 
from certain used crankcase oils. The chemical treatment, used in conjunction 
with adsorption, either prior to or after contacting with active clay, sometimes 
produces better results than can be obtained with successive applications of 
either procedure alone, thus indicating that the new and the conventional 
processes may be partially complimentary. 


LABORATORY—EXPERIMENTAL AND RESULTS 


Of the various possible selective agents, metallic sodium appeared to be the 
most effective and convenient, and the following procedure was eventually 
adopted. The oil was heated to about 120°C., sodium that had been emulsified 
previously in an additional portion of oil by mechanical agitation, at a tem- 
perature slightly above its melting point, was then added in quantities of 0.2 
to 5% by weight, and the dispersion heated to 150—200°C. Air, dispersed in 
the form of fine bubbles by a porous glass disk, was then passed through the 
emulsion for periods of time ranging from 15 min. to several hours. During this 
treatment any light ends remaining in the oil may distil from the batch and 

1Manuscript received June 13, 1955 
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some sludge will form. The amount of sludge produced varies from sample to 
sample, the smallest quantities being formed in new oils and a much larger 
amount in the case of used oils. 

The sodium is consumed in the process apparently owing to the formation 
of a complex with some of the constituents of the oil. The rate of bleaching is 
materially affected by the state of subdivision of the sodium, the rate increasing 
with decreasing particle size. The oxidation products are dispersed in the oil 
in a rather stable fashion and their removal is facilitated if coagulation can be 
induced. The most suitable method found for doing this was the dispersion 
into the system of from 1 to 25% of an oil-insoluble polar liquid such as water, 
glycerol, or alcohol. These liquids apparently wet the dispersed particles of 
sludge preferentially, causing them to collect at the oil-water interfaces. Either 
gravity settling or centrifuging then gave a satisfactory separation. These 
operations were usually performed at about 70°C. to take advantage of the 
reduced oil viscosity to promote more rapid sedimentation of the sludge. 

The oils from which the sludge has been separated appear to have a much 
improved resistance to further oxidation when compared to the unbleached 
oil. This may be due, in part, to a small amount of residual sodium complex 
that remains in solution in the oil. The presence of residual complex is indicated 
by the small increase in ash value found in the treated oil. If the oxidation is 
carried out for an extended period of time the color of the oil is found to pass 
through a minimum after which it increases again. The time required to reach 
the minimum in color is temperature dependent. With treatment temperatures 
up to 200°C., no increase in color was observed, however, as long as free sodium 
remained in the mixture. The rate of bleaching is markedly increased at higher 
temperatures and an increased amount of bleaching is usually obtained by 
using larger amounts of sodium and longer times. The treatment caused no 
substantial change in viscosity of the oils except where the optimum time was 
exceeded by a large margin. While the effect on viscosity index was usually 
small, one sample of oil with viscosity index 40 improved 10 points during the 
bleaching process. A surface active material which behaved like soap when 
titrated with acid could be recovered from the sludge. 

If the detergent properties of the oils are assessed in the manner suggested 
by Garner (2) it is found that dispersions of standard micronex (Binney and 
Smith Co. New York) are stable after some months when prepared with the 
oxidatively bleached new oil, whereas coagulation and sedimentation of the 
carbon black occur within a few days when the same oil either untreated or 
bleached with clay is used. 

Table I gives typical results of some of the laboratory tests where about 
200 cc. portions of an unbleached test oil SAE #30 were treated as described 
above, using 1% of metallic sodium. The colors are reported as percentage 
light transmission as determined by a Lumitron colorimeter (Photovolt Corp. 
New York). The instrument was adjusted to give a transmission of 100% with 
a white medicinal oil. The colors of several commercial lubricating oils, ob- 
tained in sealed cans on the open market, are shown for comparison as A-G. 
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TABLE I 


THE EFFECT OF TEMPERATURE AND TIME ON THE OXIDATIVE BLEACHING OF A 
LUBRICATING OIL 


Treatment % Transmission 


Time of oxi- 
Temp., °C. dation (min.) Red filter* Green filter* 


— None 66 28 
130 30 89 59 
130 120 94 72 
130 240 99 86 
130 720 99 87 
160 15 89 62 
160 30 99 87 
180 15 97 84 
200 5 94 72 


A 54 16 
B 74 41 
t 71 30 
D — 50 13 
E - 64 26 
F 68 34 
G 76 36 





*Red: ) = 6200A.__ 
Green: } = 5200 A. 


CONDENSER 
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4 air 
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Y 
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y 
yj 
Y 
Y 
Y 
Y 
y 
Y 
Y 
GY 


\\\\\ lato | ROTAMETER 


PUMP 


SETTLING SAMPLE 
TANK 


Fic. 1. Pilot unit. 
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PILOT PLANT—EXPERIMENTAL AND RESULTS 
In order to secure data on the behavior of larger quantities of oil, a pilot unit 
was set up that would handle four gallon batches. This is shown in Fig. 1 and 
Fig. 2. The reactor was a 6 in. diameter flanged glass vertical column, 6 ft. long. 
At the base was a metal cone in which was inserted a 2 in. diameter sintered 


* 
. ? 


Fic. 2. Pilot unit. 


stainless steel diffuser plate. The desired temperature was obtained by cir- 
culating the oil through a heated pipe. Possible overheating caused by any 
exothermicity in the reaction could be controlled by means of a cooler in the 
oil circulating line. This cooler was never required in actual practice, however. 
Light ends or water vapor that were distilled from the reacting mixture during 
a run were recovered by a condenser on the air outlet line. 

The experimental procedure was essentially the same as described for the 
laboratory work. A weighed portion of sodium emulsified in oil at 125°C. by 
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means of a Waring Blendor was dispersed in a known quantity of oil at a similar 
temperature in the reactor. The mixture was heated to the desired reaction 
temperature while air at the rate of about one cubic foot per minute was 
bubbled through the column. Oil samples were withdrawn periodically, the 
sludge coagulated with water, and the oil recovered by centrifuging at 70°C. 
In order to establish the optimum time of treatment, the reaction was fre- 
quently taken beyond the point of minimum oil color. Overtreatment normally 
had an adverse effect on the rate of settling and on the sharpness of separation. 
This effect is reflected in some of the recovery data presented, and is illustrated 
in Table II where no gravitational settling was possible in runs Nos. 1 and 6. 


TABLE II 
BLEACHING RESULTS ON NEW TEST OIL 











I ‘ 4 5 

Oil weight, Ib. 32. 32.85 32.12 30.76 
Feed 

Sodium, % 0.5 : 1.0 1.0 2.0 

Reaction temp., °C. 200 200 200 200 

Reaction time, min. 25 60 10 180 
Yield clear oil, % 

Settled 95.2 96.5 

Centrifuged 2.9 2.2 1.9 

Total 98.1 98.4 
Light transmission 

Optimum red, % 86 88 92 

Optimum green, % 43 53 69 

Time to optimum, min. 20 10 5 


Light transmission ‘Gon 40 


of original oil, % Green 10 





In these runs the time to optimum color was greatly exceeded. The optimum 
was usually indicated approximately by the color of the oil-water emulsion in 
the sludge coagulation step and by the appearance of a small amount of water 
in the condensate. 

About 15% of water gently emulsified into the oil after it was removed from 
the reactor coagulated the sludge sufficiently to give satisfactory gravity 
settling. The clear oil was recovered by decantation. Oil that was missed by 
this method of separation or that was retained in the sludge was estimated by 
centrifuging an aliquot of the material remaining in the settler after it had been 
homogenized. 

The results from the pilot unit are presented in Tables II and III. In Table 
II the oil treated was a new test oil normally used for evaluating decolorizing 
adsorbents. This oil contained no constituents volatile under the conditions of 
the procedure and it was used as received. In Table III the material described 
was a used crankcase oil obtained in 45 gal. lots from a local service station on 
two separate occasions. This oil was allowed to settle until it was reasonably 
free of sludge before treatment, and on one occasion the four gallon sample 
was washed with a 5% solution of NaOH, prior to use. About 20% of light 
ends contained in this oil, probably as the result of crankcase dilution, was 
recovered during the treatment. 
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A series of runs designed to show the effect of air flow on the bleach indicated 
that the rate of decolorizing increased as the air flow increased and as the 
pressure at which the reaction was carried out was increased. Treatment of 
samples from two other service stations (Feed Types Nos. 3 and 4 in Table III) 
indicated that the behavior reported is probably fairly representative. Re- 
moval of the volatiles prior to treatment with air and sodium had virtually 
no effect on the subsequent bleaching process. This is indicated in run No. 9 
in Table III. 

DISCUSSION 


A comparison of Table I with Tables II and III shows that colors normally 
found in commercial oils may be obtained quite easily from used crankcase 
oils, while from the new unbleached oil tested, considerably lighter oil could be 
produced with similar treatments. Recoveries from the treated new oil were 
excellent and only a small amount of sludge was formed. In the case of the used 
oils the sludge formation and sodium requirements were higher than for new 
oil. Considerable reduction in the sodium requirement was effected, however, 
with a caustic prewash and the smaller quantity of sludge formed is reflected 
in the amount of oil recovered by gravity settling. The most convenient tem- 
perature at which to work appears to be 200°C. 

For the new test oil treated with 1% of sodium at 200°C. the optimum contact 
time was in the vicinity of 20 min. For a single treatment 1% of sodium appears 
to be the optimum quantity. These points are illustrated in Figs. 3 and 4. 
Longer periods of treatment seem desirable in the case of used oils. In this case 
too an alkali prewash appears to be useful. 


o 
°o 


LIGHT TRANSMISSION % 
> 


i) 
° 


0 
0 20 40 60 
HEAT UP—>\+—TIME AT 200°C,,MINUTES 


Fic. 3. Effect of time on the decolorizing of test oil at 200°C. using 1% sodium. Curve 1, 
using red filter; curve 2, using green filter. 
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LIGHT TRANSMISSION % 


1.0 1.5 
SODIUM CONCENTRATION % 


Fic. 4. Effect of sodium concentration on the optimum color obtained from test oil at 200°C. 
Curve 1, using red filter; curve 2, using green filter. 


The chemical reactions involved in the process remain obscure. Hydrocarbon 
oxidations are generally considered to proceed by a free radical mechanism (1) 


and it would appear that in the case of the colorless species the sodium retards 
the reaction in the early stages, possibly by reacting with hydroperoxides. 
In any case the colored molecules are more susceptible to attack during the 
oxidation in the presence of sodium. 
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STATISTICAL FACTORS IN TOWER PACKING! 


By A. I. JoHNSON AND J. M. L. BEECKMANS 


ABSTRACT 


This investigation was made to determine to what extent flooding rates in 
packed extraction towers could be affected by random factors in the arrangement 
of the packing. It was found that such an effect does exist, and that the data 
could be correlated. The forms of the equations were predicted by an analysis 
based on certain simplifying assumptions. It was concluded that no exact 
prediction of flooding rates is possible, but the accuracy improves with the 
height of the packing. 


INTRODUCTION 


A knowledge of flooding rates is required for the economic design of packed 
liquid—liquid extraction columns, and a number of correlations have appeared 
in the literature. Although many refinements have been introduced, a fairly 
considerable spread of the data about the correlation equation could not be 
avoided. All these equations show that the flooding rate is independent of the 
height of the column, and that the throughput for a given ratio of disperse 
phase rate to continuous phase rate is strictly proportional to the cross- 
sectional area of the column. Nevertheless, doubts have been expressed as to 
the accuracy of the equations because of the possibility of random bottle- 
necking in the packing. Dell and Pratt (1) demonstrated that the flooding 
rate could be appreciably altered by the method used in dumping the packing. 


If air is bubbled through the column whilst it is being loaded, the flooding 
rates are reduced. 

A straight line results if the square roots of the respective phase rates are 
plotted against each other. Dell and Pratt (1) report that when a set of runs 


was repeated after redumping the packing, the straight line obtained was 
parallel to that plotted from the first set of runs, but slightly displaced. 
This is a very significant result because it can only mean that a random 
variable is present, assuming of course that the physical properties of the 
phases have remained constant. 


THEORETICAL DISCUSSION 


A randomly packed column is a highly complex geometrical system which 
defies rigorous mathematical description. The following treatment is based 
on an idealized model. It is by no means the only method of attack, but it 
has the advantage of simplicity combined with a fair degree of approximation 
to reality. 

Consider a plane section through a column, at right angles to its axis. 
Part of the area of the section will be taken by the packing, the remainder 
being available for the flow of the phases. Corresponding to this free area 

1Manuscript received June 27, 1956. 


Contribution from the Department of Chemical Engineering, University of Toronto, Toronto, 
Ontario. 
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there will be a maximum throughput which we will assume to be proportional 
to it, and independent of conditions above or below it. If the continuous 
phase rate is constant, there will be a maximum disperse phase rate which 
we will designate by v. This is the observed flooding rate in a column of 
infinitely small height which has in theory been sectioned out of a column of 
normal height. All other variables being constant, v will be a function of the 
free area of that section only. Hence v will vary statistically along the length of 
the column, and we may postulate that a distribution function of the v’s 
exists. However, owing to the finite size of the packing, v cannot vary infinitely 
quickly with the height, which we will call z. We will assume that dv/dz = 0 
over a small but definite element of height a. This does not occur in practice, 
but it will serve as an approximation to actual conditions. Fig. 1 shows the 
nature of the approximation. 


Actual distribution 


Assumed distribution 


a 
z 


Fic. 1. Typical section velocity distribution in a column. 


Let the column be divided into ” such segments of width a. The distribution 
function of the v’s will be denoted by p(v). The cumulative distribution 
function is 


P(v) = ‘a p(v) dv. 


For any one segment, P(0O < V < v) = P(v), V representing the value of 
an observation and v the value in the function. 

Since P(O < V < ~) = 1, it follows that P(v< V < ~) = 1—P(v). 
We now make the assumption that all the v’s are stochastically independent. 
In other words the probability of getting a particular v in any segment is not 
influenced by the v’s in the adjacent segments. Under these conditions the 
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probability that V should be greater than a specified value of v in m segments 
of a column simultaneously is given by the following equation: P(v < V < o) 
= [1—P(v)]". This is essentially a restriction on the minimum value of V 
in a column of height ma. Thus the right-hand side gives the cumulative 
distribution function of the minimum flooding velocities as a function of the 
number of segments and of the distribution function of the flooding velocities 
in the individual segments. If it is assumed that the flooding velocity in a 
column is equal to the flooding velocity of that section which has the minimum 
area somewhere in the column, then the above function must also be the 
cumulative distribution function of the observed flooding rates. 

Differentiation of the above equation gives the distribution function of the 
flooding rates 


[1] p*(v) = n[{1—P(v)]"""p(2). 


The next step is to determine whether the elemental height a is an inherent 
property of the packing, or whether it is possible to derive an equation be- 
tween it and the parameters of the distribution function of the elemental 
segments, in such a way that an equation may be found between 2, p*(v), 
and v. If such an equation exists, then p*(v) must be independent of a, for 
fixed values of v and z, since @ may then be eliminated. If z/a is substituted 
for m in equation [1] and the equation is differentiated partially with respect 
to a, there results 


ap(v) 
0a 


ap"(v) ee ee | ee (2a) Jay Bry ,)] (2-4) /a 
©) - 4 pei Pe)" +21-Pe)] 


2 


—*; p(v)[1— P(e)" log-[1— P(e) ] 


2(z—a) 


+2259 pip) EO) 


* [1 —P(v)|° /a@ b 


The terms in z do not drop out when this expression is equated to zero; 
hence there is no relation between a and the functions p(v) and P(v) which is 
independent of the height. 

The fundamental equation [1] is not suitable for analysis of the data in its 
present form, because nothing is known about the form of p(v). It would 
doubtless be possible to calculate the p(v) curve from values of p*(v) if suffi- 
cient data were available. However, the mathematical difficulties are con- 
siderable, principally because equation [1] is a differential equation. An exact 
solution would involve the measurement of slopes, and when one remembers 
that one is dealing with an empirical distribution function, one soon realizes 
that a solution with a relatively small number of readings with which to con- 
struct a distribution curve is out of the question. 

A better approach consists in assuming a form for p(v), with a number of 
parameters which can then be evaluated from the data. The form of p*(v) 
will depend chiefly on the form of p(v) in the small range of v. This is because 
the larger m, the greater the chances that a small value of V will occur in a 
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given column. Since one is measuring the smallest V which occurs, it is evident 
that the form of the curve p(v) at its lower end will almost entirely determine 
the whole of the p*(v) curve. This allows one to select a simple function to 
represent p(v), since the exact shape or form is important in only a relatively 
small range of the curve. The function selected is: 


[2] p(v) = C(v—v)?. 
Here v is included because p(v) must be zero at some finite value of v, since 
it is not possible to block up completely a smooth circular section with irregu- 
larly shaped areas of finite size. 
If this expression is integrated and the boundary condition that P(v) = 0 
when v = v is inserted, the following equation results: 
Cc 041 


[3] P(v) = é+1 (v—v) ° 


Then p*(v) can be obtained by substituting [2] and [3] in [1]: 


[4] p*(v) = cu 1S (o—n0)"** Fone)" 


The next step is to determine the mean of this distribution. As it stands, 
however, this is impossible, because it is only valid in the low range of v. 
Strictly the integral of the above distribution functions should tend to 1 as 
v tends to infinity. This is not so as can be seen by inspection, since the inte- 
grands themselves tend to infinity. Most distribution functions which are not 
too pronouncedly skew have a mode which is close to the mean. The mode is 
defined as that value of the independent variable for which the value of the 
function is a maximum. Furthermore the present distribution tends to get 
thinner and higher as » gets larger. This will tend to make the mode and the 
mean coincide. Fig. 2 shows the trend of p*(v) as m gets larger. 


Wit | | 


/ 


/ 


Fic. 2. Distribution of the minimum velocities. @ = 1, C = 1. 
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The mode is given by the equation dp*(v)/dv = 0, 


dp*(v)/dv = Cn(n—1) (0-09) 1 5 (o—n0)° |" (—C)(0—n" 


i. 0 6-1 
+1 +1 (v—v) O(v—v) = 0. 
The solution 1—[C/(@+1)](v—vo)*t! = 0 is trivial, since it is equivalent to 
P(v) = 1. The remaining solution is: 


bs Bice, ee 
(v—v0) ~ C{[n—1+6/(0+1)] 


If one assumes that ” > 1/(@+1), this reduces to 


0 1/(8+1) 
[5] v—09 = es . 


Or, calling 1/(6+1) = s, 
ee iw 
v0 Cz . 


Thus the mean and the height of the packing are related by an equation of 
the form: 


[6] E—vo = K/z’. 


The same difficulties are encountered in trying to evaluate the variance of 
the distribution as were encountered in estimating the mean. It is not possible 
to attack the problem directly because the appropriate integrals do net tend 
to finite values as v tends to infinity. Again this is because p(v), in the form 
assumed here, is only valid in its lower range. A qualitative estimate of the 
variance may be made by noticing that since the integral of any distribution 
function over its entire range must always be equal to 1, tall distribution 
functions will always be thin, and low ones will tend to be broad. For the 
normal distribution function there is a nice relationship between the value 
of the maximum ordinate and the standard deviation, viz., maximum 
p(v) = 1/V 2x o. 

A relationship of this form will be assumed here. First the value of the 
maximum p*(v) must be calculated. This is done by substituting equation 
[5] into equation [1]. 


6 6 /(6+1) 
Thus = Cn ’ 


If it is assumed that 7 is large, the following approximation may be made: 
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where e is the base of the natural logarithms. Substituting and simplifying, 
one obtains 


e 6 /(8+1) 
[7] p*(v)m = cul #) ‘ 


If one makes the assumption that this value is inversely proportional to the 
variance, as in the normal distribution, then the standard deviation and the 
height must be related by an equation of the form 







[8] o = K's’, 






The exponent ¢ is related to the exponent s of equation [6] by the following 
equations: 







0 aa 1 
i ~(:-;4) = 941 2S = FT 





Hence 


[9] t= -—s. 








A small correction may be brought in at this stage. The observed standard 
deviations are due to two causes, one being that of the random elements in 
the packing, the other being the error in making the measurements. The 
theory of probability states that the standard deviation of a variable which is 
the sum of two independent variables is related to the standard deviation of 
these variables by the equation: 


[10] c= o;°+0-2?. 










If one postulates a distribution function of the experimental errors, with 
standard deviation ¢,, then the following equation gives the relation between 
the observed standard deviation § and that due to random factors in the 
packing: §2 = o?+o,?. The relation between § and z then assumes the form 
8’—o¢,2 = K’z-’. This concludes the derivation of the correlating equations. 







APPARATUS AND PROCEDURE 

The arrangement of the apparatus is shown schematically in Fig. 3. The 
column had an internal diameter of 6 in. All conduits were of 13 in. copper 
tubing. The method of supporting the packing and the connections between 
the column and the end-pieces are shown in Fig. 4. The packing consisted of 
2 in. ceramic Raschig rings. 

The solvent used was a mixture of carbon tetrachloride and gasoline, the 
proportions being such that the mixture had a density of 1.340. 

In operation the procedure was to drop in packing to a certain height, 
determine the flooding rate, then add more packing, etc. The flooding rate 
was taken to be that rate which caused a layer of solvent about 6 in. wide 
above the packing neither to increase nor to decrease with time. Flooding 
was observed at packing heights of 3 in., 6 in., 1 ft., 1 ft. 6 in., 2 ft., 2 ft. 6 in., 
3 ft., and 4 ft. At the end of a series of runs the column was repacked and so on. 
A total of 19 series were made. A close check was held on the density by means 
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Fic. 3. Schematic flow diagram. 


Fic. 4. Section through base of column. 
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of a hydrometer. The permissible limits of variation were 0.010 about a mean 
of 1.340. 

It should be noticed that the readings in a series (Table I) are not indepen- 
dent of each other, since the column is not entirely repacked at each reading. 
However, in a fairly large sample such as this, this could not be of importance. 

Unfortunately, parts of the column corroded during the experiments, 
and corrosion products contaminated the phases. In particular, it is thought 
that the interfacial tension was affected. As a result, the readings at any 
given value of the height tended to rise as the investigation proceeded. 

Fig. 6 shows that the data can be correlated by an equation of the type 
predicted, equation [6]. The value of the exponent s is 1.284, and v has a value 
16.5 lb./min. 


TABLE I 


FLOODING RATES, LB./MIN. 
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TABLE II 


CALCULATED VALUES OF THE MEAN AND STANDARD DEVIATION 


=x =x? Ex a 


534.1 16266 .8 15846 .7 24.64 

448.2 ; 11429 .1 11160.2 15.81 

390.6 : 8120.3 8026.8 5.022 
356.8 6772.9 6697 .1 4.033 
345.5 6339 . 4 6281.2 3.150 
335.8 5991.1 5933 .6 3.127 
330.5 5786.8 5747.4 2.100 
326 .6 5642.5 5617.5 1.583 
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The same can be said of the correlation of the variances, shown in Fig. 7, 
corresponding to equation [8]. The parameters are o,? = 0.3 and s = 0.532. 

The agreement between the two values of s is not good, but this is to be 
expected in view of the approximations made and the experimental error. 
Fig. 5 shows plots of the flooding rate versus chronological sequence of the 
readings for packing heights of 3 in. and 4 ft. The trend is obvious, even 
allowing for legitimate fluctuations. It is interesting to notice that the drift is 
much more pronounced at low heights of packing. Although the equations 
developed are now no longer strictly applicable, the qualitative predictions 
remain true nevertheless. 


CONCLUSIONS 


The agreement between theory and experiment is as good as could be 
expected, and it is reasonable to conclude that a mechanism such as that 
postulated is operative. The percentage deviation at 4 ft. of packing, 7.3%, 
is quite small considering that it is inflated owing to the experimental error. 

It can therefore be concluded that statistical factors are not very important 
in affecting the flooding rate in a column of this kind, using 4 in. Raschig 
rings, with a packing depth of 4 ft. or more. This conclusion cannot be extended 
to other cases without further experimental work. It is expected that the ratio 
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of packing size to column diameter would be an important parameter affecting 
the statistical effect, as well as the ratio of packing size to depth of packing, 
and the type of packing used. It is recommended that no experimental work 
be done on columns for which the above ratios are greater than in the present 
case. 


NOTATION 
a Height of an elemental segment, ft. 
C A constant defined by equation [2]. 
s Exponent defined by equation [5]. 
§ Standard deviation of the flooding rates at constant height of packing, 
lb./min. 
t Exponent defined by equation [8]. 
v Theoretical flooding rate in the elemental segment, Ib./min. 
vo Smallest possible flooding rate, lb./min. 
V Observed flooding rate, lb./min. 
6 Exponent defined by equation [2]. 
o Theoretical standard deviation, Ib./min. 
& Mean flow rate at a given height, lb./min. 
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THE OXIDATION, DECOMPOSITION, IGNITION, AND 
DETONATION OF FUEL VAPORS AND GASES 


XXVII. THE HYDROGEN ENGINE! 


By R. O. Kinc? AND MoGeEens RAND? 


ABSTRACT 


Experiments made with hydrogen as the fuel for an Otto cycle engine show 
that optimum performance in terms of power and economy is limited only by 
preignition arising from the temperatures attained by the ceramic core of the spark 
plug or by the exhaust valve, or from nuclei formed by pyrolysis of lubricating 
oil vapor. On reduction of the magnitude of these igniting ‘effects it became 
possible to operate the engine with speeds rising to 1800 r.p.m. and the com- 
pression ratio fixed at 12:1. A value of 146 p.s.i. for I.M.E.P. with an indicated 
thermal efficiency of 37% was obtained when the maximum power mixture 
strength was used. A value of 90 p.s.i. was obtained for I.M.E.P. with an 
indicated thermal efficiency of 47% with a mixture 55% weak. These values 
exceed, by a wide margin, any recorded in the literature of the subject. 


INTRODUCTION 


Results in conformity with the nuclear theory of ignition, as stated in 
Part IV (7), were obtained when hydrogen was used as the fuel for experiments 
with a C.F.R. knock rating engine (Part V (12)). These experiments were 
carried out at the standard knock rating speed of 900 r.p.m. and with com- 
pression ratios rising to 10:1. Ignition otherwise than by the electric spark did 
not occur unless finely divided carbon had been allowed to accumulate in the 
combustion chamber. It was found later that the carbon was a product of the 
pyrolysis of lubricating oil vapor. A similar igniting effect occurred when 
charcoal dust was added to the hydrogen—air mixture before it entered the 
combustion chamber and when in subsequent experiments the laboratory 
atmosphere was impregnated to an unusual degree with concrete dust (11). 
This was of especial interest in that the nuclei were not combustible. A signifi- 
cant feature of the experiments was the occurrence of some degree of pre- 
ignition on increase of the temperature of the jacket coolant from 140° to 
212°F. This effect according to the nuclear theory of ignition would be due toa 
surface in the combustion chamber having attained an igniting temperature 
or to nuclear ignition arising from an increase in the rate of vaporization of 
lubricating oil. The surfaces that would be expected to attain an igniting 
temperature were those of the exhaust valve, and the ceramic core of the spark 
plug or of its central electrode. The most effective way of raising the tempera- 
ture of these indirectly cooled surfaces would have been to increase the heat 
load by increasing engine speed but the unbalanced C.F.R. knock rating engine 


1 Manuscript received May 31, 1955. 
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ment of Mechanical Engineering, University of Toronto, Toronto, Ontario. 
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was not suitable for operation at speeds higher than 900 r.p.m. The engine 
experiments were therefore discontinued until a C.F.R.-F4 engine, balanced 
for high speed operation, became available for the new series of experiments 
described in this Part. 

Experiments with combustion tubes of materials similar to those in the 
combustion chamber of the engine were however initiated in order to determine 
the wall temperature required for the ignition of hydrogen—air mixtures passing 
through them and are described in Part XXVI (10). 


EXPERIMENTAL ARRANGEMENTS AND METHODS, PRESENTATION OF 
RESULTS 
The Engine 

The C.F.R.-F4 engine, bore 3} in. stroke 43 in., is fitted with balancing 
pistons in the crank case, as required for high speed running. The cast iron 
piston is fitted with an oil scraper ring and three pressure rings. The engine 
was direct connected to a d-c. generator arranged as a swinging field dyna- 
mometer for the measurement of B.H.P. directly and of I.H.P. by the motor- 
ing method. Provision was made for the measurement of compression ratios 
varying from 4 to 20:1. 

Camshafts which provided suitable valve timing for speed ranges of 900 
to 1200 r.p.m. and 1200 to 1500 r.p.m. were available. They will be described 
respectively as low and high speed camshafts. The valve timing and lift are 
given in the appendix to this Part. 

Cylinders with holes for three spark plugs in the wall and one in the head 
were available. One was fitted with a sodium cooled exhaust valve, and others 
with the C.F.R. standard uncooled valve of silchrome steel. The types of the 
cooled and uncooled valves are illustrated by Fig. 15 of the appendix. A plain 
or a shrouded inlet valve could be used in either cylinder. 


Exhaust Valve Surface State 


The surface state of the face of the exhaust valve varies with time of use 
and this was found to have an appreciable effect on the temperature required 
for the ignition of other than weak hydrogen-air mixtures. Exhaust valves 
with two varieties of surface state were used, namely clean and aged. 

A valve cleaned with abrasive material and then washed with alcohol before 
a set of experiments was begun is described as clean although it ceases to re- 
main in that condition to an extent depending on the time required to complete 
a set. A running time of from three to four hours is generally required. 

An aged valve is one that has been used for many sets of experiments and 
cleaned only by washing with alcohol, before a set of experiments is begun. 


The Ignition System 


The designs of spark plugs used for the experiments are illustrated by 
Fig. 14 of the appendix. Champion C.F.R.-8 was until recently the standard 
spark plug for use with the A.S.T.M.-C.F.R. method of knock rating. It is a 
“hot” variety, near the top of the commercial heat scale. Spark plugs K.L.G. 
100 and Champion R 2-S were used as cool varieties. 
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The coil ignition supplied as standard had been found to be ineffective with 
hydrogen—air mixtures at a C.R. of 10.0, unless the spark gap were reduced to 
0.012 in. A special Auto-Lite ignition system was used for experiments made 
subsequently with hydrogen and other gases and it was then possible to in- 
crease the spark gap to 0.025 in. if necessary. The system is described in the 
appendix to Part XIII (9). 


Temperature Control 


It was found in the course of preliminary experiments that as the heat load 
carried by the engine was raised, by increasing speed for example, it became 
necessary to control the temperatures of the air and hydrogen supplied to the 
engine if repeatable results were to be obtained. The temperature of the air 
drawn from the laboratory atmosphere was subject to variations of approxi- 
mately 40°F. The method used for maintaining the supply of air and hydrogen 
at a nearly constant temperature of 65°F. is described in the appendix to this 
Part. 

Water was used as the jacket coolant and temperatures lower than 212°F. 
were controlled automatically by the thermostatic method described in Part 
XIII (9, pp. 152-153). The flow of water through the cylinder jacket could 
be controlled by this method to ensure that the temperature at the outlet 
was never more than 5°F. higher than that at the inlet. 

The lubricating oil in the crank case sump was maintained at 140°F. by an 
exterior electric heater with manual control. Oil pressures maintained at 
speeds of 900, 1200, 1500, and 1800 r.p.m. were 18, 21, 23, and 25 p.s.i. respec- 
tively. 


The Hydrogen Supply 

Electrolytic hydrogen as supplied by the Canadian Liquid Air Company, 
Toronto, compressed to 2000 p.s.i. in steel cylinders was used for the experi- 
ments. The hydrogen was reputed to be 99.95% pure and as compressed was 
nearly dry, the dew point being below minus 100°F. The lower calorific value 
as used for calculations of thermal efficiency was taken as 274.5 B.T.U. per 
cu. ft., at 60°F. and 30 in. Hg. The method used for the measurement of the 
rate of hydrogen supply to the engine is described in the appendix to this 
Part. 


The Lubricating Oil 


A separate series of experiments with lubricating oils indicated that the 
least amount of residue remained in the combustion chamber after Ucon had 
been used. This is a synthetic lubricant made by the Carbide and Chemicals 
Corporation with a base composed of polyalkylene glycols and derivatives. 
The grade selected for the experiments is described as L.G. 300 X, which 
indicates that it is water insoluble, has a viscosity of 300 Saybolt universal 
seconds, and is oxidation inhibited but without additives. According to 
A.S.T.M. tests, the viscosity index is 142, ash less than 0.01%, carbon residue 
also less than 0.01%, and water less than 0.25%. The viscosity is close to 
S.A.E. 30 at 210°F. and the same as S.A.E. 20 at 130°F. 
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Methods of Experiment 


The engine was always prepared for a particular set of experiments by being 
dismantled for a top overhaul. The cylinder, piston, and rings were then 
washed with acetone and alcohol. The valves were inspected and hand lapped 
into their seats if signs of leakage were detected. On conclusion of the set it 
was verified, by measurement of compression pressure before the engine had 
cooled, that piston sealing had not deteriorated. 

The lubricating oil and the jacket cooling water were heated electrically to 
working temperatures before a set of experiments was begun. The engine was 
then started on a hydrogen—air mixture 50% weak and run for at least 30 min. 
The mixture was then reduced to the leanest on which the speed selected for the 
first experiment could be maintained. This was always the lowest of the 
speeds to be used for a set of experiments. The engine was then run for 15 
min. before a determination of the optimum spark timing. 

Hydrogen burns with high velocity in mixtures with air approaching and 
exceeding the correct value and relatively slowly when in small concentration 
in air, and the ignition timing required for maximum power output when a 
particular mixture strength and engine speed was used can be determined with 
accuracy, solely by a series of lower measurements. Brake horsepower was 
measured accordingly for a series of ignition timings giving rising and falling 
values for power output and the optimum timing was then obtained from a 
plot of results, as described in Part XII (8, pp. 488-439). 


Presentation of Experimental Results 


Experiments were carried out at engine speeds rising to 1800 r.p.m. and with 
compression ratios of 10 and 12:1. Comprehensive performance data are given 
by graphs plotted on a mixture strength base as used for experiments with 
liquid fuels. A ‘‘correct’’ mixture is 100% stoichiometric. The percentage 
stoichiometric diminishes or increases as the mixture is weakened or enriched. 
Thus for example a mixture 20% weak is 80% of stoichiometric and one 20% 
rich is 120% of stoichiometric. 

The terms Compression Ratio, Mixture Strength, Indicated Thermal 
Efficiency, and Indicated Mean Effective Pressure are, when convenient, 
abbreviated to C.R., M.S., I.T.E., and I.M.E.P. respectively. Values of C.R. 
are given as single numbers. A stoichiometric mixture is described as ‘‘correct”’ 
in the text; the abbreviation C.C.M. is used on the graphs. It is to be under- 
stood that numerical values given on graphs for I1.M.E.P. and Engine Speed 
are in pounds per square inch (p.s.i.) and revolutions per minute (r.p.m.) 
respectively. 


I. THE IGNITING EFFECT OF THE CERAMIC CORE OF A SPARK PLUG 
Preliminary Experiments 
The operating conditions were similar to those in which some degree of 


preignition occurred in the course of the experiments described in Part V 
(12) when a correct hydrogen—air mixture was used at a C.R. of 10, with the 
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jacket coolant at 212°F. and the mixture supplied at room temperature. A 
C.F.R.-8 spark plug which is a “‘hot’”’ variety was used, the exhaust valve was 
of silchrome steel, uncooled, and the inlet valve was of the shrouded type 
which provided a swirl in the entering mixture. Thus the conditions were 
essentially those used for the knock rating of fuels by the C.F.R.-A.S.T.M. 
motor method. 

The experiments were carried out with the C.R. fixed at 10, with engine 
speeds of 600, 900, 1200, and 1500 r.p.m. while mixtures were varied from 
very weak to extremely rich. The results obtained for the relation between 
I.H.P. and M.S. are given by the graphs of Fig. 1. 


| 
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Fic. 1. Mixture strengths over which premature ignition occurs when engines speeds 
exceed 900 r.p.m. and a C.F.R.-8 spark plug, a hot variety, is used in the standard position. 


It will be seen by reference to the graphs that at 600 r.p.m. a maximum 
value of 2.6 was obtained for I.H.P. when the mixture was 10% rich. There 
was no preignition. When however the heat load was increased by raising the 
speed to 900 r.p.m. and I.H.P. attained a maximum of 4.3, there was some 
degree of preignition but not enough to interrupt the experiment. This result 
was similar to that obtained in similar conditions as described in Part V (12). 
Then when the heat load was further increased by raising the speed to 1200 
and then to 1500 r.p.m. the range of M.S. over which running was impossible, 
because of backfiring into the induction system, increased as shown by the 
relevant graphs. Thus at 1500 r.p.m. backfiring occurred over the mixture 
range 40% weak to 140% rich. 


Ignition by a Hot Spark Plug without a Central Electrode 


Preignition is frequently attributed to the central electrode of the spark 
plug becoming overheated. Pye (15, p. 91), for example, mentions only the 
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central electrode and scraps of incandescent carbon; Hundere and Bert (6, 
p. 553), fragments detached from the ceramic core. The experiments described 
in Part XXVI (10) indicated that preignition would preferably be initiated by 
the surface of the ceramic core of the spark plug even if its temperature were 
not higher than that of other surfaces in the combustion chamber. 
Experiments were carried out accordingly in the operating conditions used 
for the preliminary experiments but with the mixture fired by a cool variety of 
spark plug, a K.L.G.-100, in the standard position while a standard but unfired 
hot variety, a C.F.R.-8, was in the No. 2 position. The central electrode of this 
plug had been cut off. The standard and the No. 2 positions are shown by 
Fig. 2. It will be seen by reference to the figure that the swirl due to the 


STANDARD POSITION 





Fic. 2. Location of spark plug holes in wall of the combustion chamber and the mixture 
swirl due to the use of the C.F.R. shrouded inlet valve. 


shrouded inlet valve directs the relatively cool entering mixture over the spark 
plug in the standard position and that it then passes over the hot exhaust 
valve before arriving at the unfired plug in the No. 2 position. 

The experimental results obtained for the relation between I.H.P. and M.S. 
with speeds of 600, 900, 1200, and 1500 r.p.m. are given by the graphs of Fig. 3. 
It was just possible to use the maximum power mixture strength without the 
occurrence of preignition at the low speed of 600 r.p.m., the corresponding 
maximum I.H.P. being 2.6. Then on raising the speed to 900 r.p.m. preignition 
accompanied by backfiring occurred with mixtures ranging from 16% weak to 
67% rich and the range of mixture strength over which running was impossible 
increased as speed and the corresponding heat load were further increased, as 
illustrated by the graphs of the figure. 
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THE HYDROGEN ENGINE 
LOW SPEED CAMSHAFT, C.R. 10. JACKET COOLANT 
212 °F. MIXTURE SUPPLY, 65 °F. CLEAN 


SILCHROME EXHAUST VALVE NOT COOLED 
SHROUDED INLET VALVE 


REGION OF PREIGNITION AND 
EXPLOSION | ww eNO UCTION 
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Fic. 3. Mixture strengths over which premature ignition occurs at engine speeds exceeding 
600 r.p.m. when an unfired hot spark plug without projecting central electrode is in No. 2 
position and ignition is by a cool spark plug in the standard position. 


II. THE LIMITATION OF POWER AND ECONOMY BY THE ONSET OF PREIGNI- 
TION DUE TO THE ATTAINMENT OF AN IGNITING TEMPERATURE BY THE 
CERAMIC CORE OF THE SPARK PLUG 

The experiments described in Section I were carried out at a C.R. of 10 and 
in the absence of an accumulation of finely divided carbon in the combustion 
chamber. The maximum power attainable was limited by preignition arising 
from the ceramic core of the C.F.R.-8 spark plug having attained an igniting 
temperature in respect of the hydrogen-air mixture. Moreover it was shown 
by the experiments with combustion tubes (Part X XVI (10)) that a correct 
mixture would be ignited by the exhaust valve if it attained a temperature of 
approximately 1100°F. It was decided accordingly to use a sodium cooled 
valve for further experiments and to maintain the jacket coolant at 100°F. 
and the entering mixture at 65°F. The shrouded inlet valve was replaced by a 
plain type although this led to an increase in volumetric efficiency and conse- 
quently to an increase in heat load. The volumetric efficiency at high speeds 
was however restricted to some extent by the continued use of the low speed 
camshaft. The object of the changes mentioned was to ensure that preignition, 
when it occurred, was due solely to the temperature attained by the ceramic 
core of the spark plug. Experiments were carried out at a C.R. of 10 and with 
speeds of 900, 1200, 1500, and 1800 r.p.m. when a hot spark plug was used and 
repeated after a change to a cool variety of spark plug. 


Hot Spark Plug Experiments 

A Champion C.F.R.-8 was used as a hot plug. The relation between I.H.P. 
and the values of M.S. that could be used before the experiment was interrup- 
ted by the occurrence of preignition and backfiring is given by the A graphs of 
Fig. 4, for the four engine speeds. Backfiring began on increasing M.S. to the 
values shown by the short vertical lines of the graphs. The mixture was then 
42% weak at speeds of 1500 and 1800 and 38% weak at 900 and 1200. A value 
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of 40% weak, taken irrespective of engine speed, would be within experimental 
accuracy. Preignition and backfiring did not occur when extremely rich mix- 
tures were used. The values then obtained for I.H.P. are not given by the 
graphs. The small ranges of M.S. usable in the conditions of the experiments 
included the weak value for which indicated thermal efficiency was a maximum. 
Relevant data are tabulated below. 


TABLE I 


MIXTURE STRENGTH FOR MAXIMUM I.T.E., 60% WEAK 
ONE HOT SPARK PLUG, LOW SPEED CAMSHAFT, C.R. 10.0 


Opt. spk. 
Speed Advance .M.E.P. 
30° 
1500 32° 
1200 32° 
900 32° 


THE HYDROGEN ENGINE C.R.10 


CLEAN SODIUM COOLED EXHAUST VALVE 
LOW SPEED CAMSHAFT, MIXTURE SUPPLY, 65°F pater COOLANT, 00°F 
T 
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¥ Fic. 4. Graphs A and B illustrate the effect of changing from a hot to a cool variety of 
spark plug; other running conditions having remained unchanged. 


Cool Spark Plug Experiments 

A K.L.G. 100 was used as a cool plug. The relation between I.H.P. and M.S. 
then obtained without otherwise changing conditions of experiment is given 
by the B graphs of Fig. 4. They show that a remarkable effect, which prevented 
premature ignition, had been obtained. Thus referring to the figure it will be 
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seen that hydrogen—air mixtures ranging from 70% weak to 80% rich were used 
at speeds ranging from 900 to 1800. 

The values obtained for I.T.E. and I.M.E.P. when the M.S. was used for 
maximum power and optimum spark timing are given, in Table II, below. 
Maximum power was always obtained when using a mixture 5% rich. The 
values for I.T.E. and the corresponding values of I.M.E.P. obtained when 
using the M.S. for maximum I.T-.E. are given in Table III. 


TABLE II 


MIXTURE STRENGTH FOR MAXIMUM POWER, 5% RICH 
ONE COOL SPARK PLUG, LOW SPEED CAMSHAFT, C.R. 10 








% Vol. Spk. 
ff. timing 


1° adv. 
t.d.c. 

1° ret’d 
1° ret’d 








TABLE III 


MIXTURE STRENGTH FOR MAXIMUM THERMAL EFFICIENCY, 55% WEAK 
ONE COOL SPARK PLUG, LOW SPEED CAMSHAFT, C.R. 10 


Speed, % Vol. Spk. 
eff. 


r.p.m. adv. LEP. I.M.E.P. 





1800 75 35° 6.1 75 
1500 85 26° 6.0 85 
1200 91 28° 4.9 87 

900 22° 3.6 85 





It is shown by the data of Tables II and III that a maximum value of 134 
was obtained for I.M.E.P. with a speed of 1200 and the somewhat lower 
value of 132 on increasing speed to 1500. A maximum of 42.4% was obtained 
for I.T.E. with a speed of 1200 and of 44% on increasing the speed to 1500. 
The I.H.P. obtained with the maximum lower mixture increased from 4.9 
to 6.0 on increasing the speed from 1200 to 1500 but an increase to 6.1 only 
was obtained on further increasing speed to 1800 and the I.M.E.P. decreased 
from 132 to 123. 

It is evident that the valve timing provided by the low speed camshaft was 
suitable for use at speeds of from 1200 to 1500 but not for a speed of 1800. 
Comprehensive performance data were obtained for the four speeds used for 
the experiments and with mixtures ranging from 66% weak to 60% rich. 
The data are given by the graphs of Fig. 5 for a speed of 1500 only. 
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Fic. 5. Performance data, mixtures 66% weak to 60% rich, obtained when ignition is by a 
cool spark plug, speed 1500 r.p.m. and C.R. 10:1. 


Ill. THE MEAN TEMPERATURES ATTAINED BY HOT AND COOL SPARK PLUG 
CORES, AN UNCOOLED EXHAUST VALVE, AND THE EXHAUST GAS AS 
AFFECTED BY MIXTURE STRENGTH 
The experiments were carried out with an engine speed of 1200 r.p.m., a 
C.R. of 10, temperatures of jacket coolant and mixture supply 100° and 65°F. 
respectively. The experimental results are given by the graphs of Fig. 6. 


The C.F.R.-8 Spark Plug 


The spark plug with a thermocouple in the central electrode near the tip 
of the ceramic core was supplied by the Champion Spark Plug Company. 
Thus the temperature at that position could be measured while the spark 
plug, placed in the standard position, was used to fire the hydrogen-—air 
mixture. The temperatures attained as M.S. was increased from 70% weak are 
given by the A graph of the figure. Preignition occurred on increasing the M.S. 
to 33% weak. A temperature of 1150°F. had then been attained at the position 
of the thermocouple. 


The Uncooled Exhaust Valve 


The valve was of silchrome steel, supplied as standard for the C.F.R.- 
A.S.T.M. motor method of knock rating. It was obtained with the stem 
drilled for a thermocouple to be fitted in a manner suitable for the measure- 
ment of the temperature of the exposed surface. The temperatures thus 
measured while ignition was provided by a cool spark plug, a K.L.G. 100, in 
the standard position, are given by the B graph of the figure. It will be seen 
that the temperature of the exhaust valve rose from 590° to a maximum of 
1150°F. as M.S. was increased from 66% weak to 5% rich. There was no 
preignition. 
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Fic. 6. Mixture strength and the corresponding temperatures of an uncooled exhaust valve, 
the exhaust gas, and hot and cool spark plugs. 


The Exhaust Gas 


The temperature was measured by using a bare thermocouple placed in the 
exhaust system approximately 3 in. from the valve. The temperatures thus 
measured, as M.S. was increased from 60% weak to 20% rich, are given by the 
C graph of the figure. They were measured at the same time as those of the 
uncooled exhaust valve. It will be noted that the temperature of the uncooled 
exhaust valve was always higher than that of the exhaust gas. The temperature 
of a sodium cooled valve in an air cooled cylinder was always lower than that 
of the exhaust gas, according to experiments by Sanders, Wilstead, and 


Mulcahy (17). 


The Cool Spark Plugs, R2-S and K.L.G. -100 


An R2-S spark plug with a thermocouple in the central electrode near the 
top of the ceramic core was placed in the standard position and used for firing 
the mixture. A K.L.G.-100 was prepared for the experiments by removing the 
central electrode and replacing it by a thermocouple cemented in at about 
1/8 in. from the tip of the ceramic core. This plug was placed in the No. 2 
position. The inlet valve was of the plain type, thus there was no swirl to 
influence the relative temperatures attained by the spark plugs. The K.L.G. 
plug in the No. 2 position was however closer to the exhaust valve than the 
champion plug in the standard positions. The temperatures, as measured while 
M.S. was increased from 70% weak to 20% rich, are given by graphs D and 
E of the figure. It will be seen that the temperatures attained by the two 
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plugs did not differ to a significant extent, allowing for their relative positions, 
and that the maximum temperature of 670°F., which was attained with a 
mixture 10% rich by the K.L.G.-100 plug, is far below that required for the 
ignition of any hydrogen-air mixture. 









IV. THE LIMITATION OF POWER AND THERMAL EFFICIENCY BY THE ONSET 
OF PREIGNITION DUE TO THE EXHAUST VALVE HAVING ATTAINED AN 
IGNITING TEMPERATURE 


























The experiments described in Sections II and III demonstrated that pre- 
ignition due to the ceramic core of the spark or to its central electrode did not 
occur when a cool type of spark plug was used. It was shown however, by the 
B graph of Fig. 6, that the temperature attained by the exhaust valve ap- 
proacned the value required for the ignition of a correct hydrogen—air mixture 
when the engine speed was 1200 r.p.m. and the C.R. 10:1. It appeared therefore 
that preignition by the exhaust valve might become a limiting factor in respect 
of power and thermal efficiency on raising the heat load. 

Experiments were therefore carried out in the operating conditions used for 
those described in Sections II and III, except that volumetric efficiency and the 
consequent heat load were increased by using a “high speed’’ camshaft. 
There was then a corresponding increase in the temperature of the exhaust 
valve and a further increase due to the valve lift being increased from 0.238 
in. to 0.312 in. and to the time the valve was off its cooling seat being increased 
from 145° to 155° of crank revolution. The values for I.H.P. obtained with 
mixtures ranging from 70% weak to 60% rich are given by the graphs of 
Fig. 7. The I.H.P. attained a maximum of 11.0 with a speed of 1800 r.p.m. as 
compared with the value of 10.4 obtained when using the low speed camshaft, 


HIGH SPEED CAMSHAFT 
ONE COOL SPARK PLUG, CR. 10 


INDICATED HORSEPOWER 
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Fic. 7. Premature ignition at speeds lower than 1500 r.p.m. as obtained when changing to 
the high speed camshaft, other running conditions remaining unchanged. 
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Fig. 4. The especially noteworthy feature of the experiments was that preigni- 
tion and backfiring into the induction system occurred when speed was 
reduced. 

An explanation was provided by the results of the experiments with steel 
combustion tubes, Part XXVI (10). It was found in these experiments that 
the igniting effect of a clean stainless steel surface was reduced if it was aged 
by exposure at high temperatures to hydrogen—air mixtures for the time 
required for the formation of an adherent layer of oxide. Moreover it was 
shown that the difference in the igniting effects increased as time of exposure 
to hydrogen—air mixtures was increased. Thus it would be expected that in the 
engine an exhaust valve with a clean surface would ignite a hydrogen-air 
mixture at a lower temperature than would be required by one with an aged 
surface and that the difference in the igniting effects would increase on 
decreasing the engine speed. 

The surface of the exhaust valve had always been cleaned with abrasive 
material for the experiments described earlier in this Part but in view of the 
results of the experiments with stainless steel combustion tubes the clean 
exhaust valve was replaced by a similar sodium cooled valve which had been 
put aside after having been used for a long period during which a coating of 
oxide had formed. The valve is described accordingly as ‘‘aged’’. 

The engine then ran at speeds ranging from 900 to 1800 r.p.m. with an 
entire absence of preignition. No alteration other than the change of exhaust 
valves had been made in the operating conditions. The values obtained for 
I.H.P. when mixtures ranging from 70% weak to 70% rich were used are given 


by the graphs of Fig. 8. The experiments were of a preliminary nature made to 
verify that preignition would be avoided by using an aged instead of an 
initially clean exhaust valve. 
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Fic. &. Graphs showing the absence of low speed premature ignition when the exhaust 
valve was aged. 
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Performance Data for Six Engine Speeds 

Further experiments were carried out to obtain complete performance data 
when the aged exhaust valve was used with speeds of 1100 and 1300 r.p.m. 
in addition to the speeds of 900, 1200, 1500, and 1800 r.p.m., while mixtures 
were varied from 70% weak to 70% rich. The data obtained when the mixtures 
required for maximum I.H.P. and maximum [.T.E. were used are tabulated 
in Tables IV and V, respectively. 


TABLE IV 


MIXTURE STRENGTH FOR MAXIMUM POWER; AGED EXHAUST VALVE 
HIGH SPEED CAMSHAFT, ONE COOL SPARK PLUG, C.R. 10 


Spk. 


% Vol. 





Speed eff. timing M.S. EELP. I.M.E.P. I.T.E.% 
1800 81 5S 2° adv. 5% rich 11.40 135 36.5 
300 «89.0stdcs «he rich «= BLOB 132033 
1100 «8718 tde, «BSE rich «6.0 128 38.5 

900 84.5 2° retard 5% rich 5.25 122 33.0 


TABLE V 


MIXTURE STRENGTH FOR MAXIMUM INDICATED THERMAL EFFICIENCY, 55% WEAK; 
AGED EXHAUST VALVE 
HIGH SPEED CAMSHAFT, ONE COOL SPARK PLUG, C.R. 10 








% Vol Spk. 
Speed efi. adv. 1.H.P. I.M.E.P. 1.T.E.% 
1800 83 24° 7.0 83 44.4 
1500 87 7° 6.1 87 43.2 
1300 84 30° 5.1 84 42.0 
1200 84 31° 4.8 84 41.5 
1100 82 30° 4.2 82 41.5 
900 80 30° 3.3 80 40.2 


It will be seen by reference to Table IV that I.H.P. was 11.4 at 1800 r.p.m., 
as compared with 10.4 H.P. obtained when the low speed camshaft was used, 
see Table II. The maximum values for I.M.E.P. and volumetric efficiency were 
obtained at 1500, indicating that the valve timing provided by the high speed 
camshaft was most suitable for that engine speed provided that preignition 
was avoided by using an aged exhaust valve. The I.M.E.P. was then 139 as 
compared with a value of 132 p.s.i. obtained when the low speed camshaft and 
a speed of 1500 was used, and 134 when the speed was 1200 r.p.m. It may be 
concluded that in so far as I.M.E.P. is concerned, the valve timing provided 
by the low speed camshaft is most suitable for an engine speed of 1200 and 
that provided by the high speed camshaft is best for an engine speed of 1500 
r.p.m. 

The maximum values for I.T.E. were always obtained when the mixtures 
were 55% weak, when either the low or the high speed camshaft was used. 
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The values obtained with the low speed camshaft, Table III, increased from 


42% at a speed of 900 to attain a maximum of 44% at 1500 r.p.m. and a lower 
value of 43.6% at 1800 r.p.m. On the other hand the values obtained when 
the high speed camshaft was used increased with fair regularity from a low of 
40.2% at 900 to a maximum of 44.4% at 1800 r.p.m., Table V. Performance 
data for mixtures ranging from 70% weak to 80% rich are given by the graphs 
of Fig. 9 for the engine speed of 1500 at which a maximum of 139 was obtained 
for I.M.E.P. 
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Fic. 9. Performance data, mixtures 70% weak to 64% rich, obtained when using the high 
speed camshaft and the aged exhaust valve. Speed 1500 r.p.m. C.R. 10:1. 


Indicator Diagrams 

The indicator diagrams, Fig. 10, were taken when a Draper-Li pressure 
pickup was used. The pressure scale is linear and maximum values in pounds 
per square inch are given by the calibration marks on the diagrams. The 
diagrams were taken when a C.R. of 10 was used and when, under the operating 
conditions used, the performance data given in Tables IV and V were obtained, 
at an engine speed of 1500 r.p.m. 

The pressure waves started by the spark are most pronounced, when with 
mixtures nearly correct the pressure rise due to combustion began immediately 
after its passage, diagrams A and B. Spark timing is shown on the diagrams 
by small rather faint blips fed into the electronic circuit. These are marked by 
arrows. Similar blips marked by short vertical lines were fed into the circuit 
at t.d.c. and at intervals of 10° for 30° before and after. 

Spark timing was always adjusted to obtain maximum power output and 
it will be noted that maximum combustion pressure then occurred at approxi- 
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Fic. 10. Pressure-time diagrams taken with a Draper-Li pressure pickup. C.R. 10:1 speed 
1500 r.p.m. 


A, M.S. 3% rich, spk. 1° retard. B, M.S. 8% weak, spk. 1° adv. 
C, M.S. 34% weak, spk. 14° adv. D, M.S. 50% weak, spk. 25° adv. 
E, M.S. 80% weak, spk. 75° adv. F, air only. 


mately 15° after t.d.c. when the mixture was nearly correct, diagrams A and 
B, but at approximately 5° after t.d.c. when mixtures 34, 50, and 80% weak 
were used, diagrams C, D, and EF. 

The power developed with a mixture 80% weak was insufficient to over- 
come friction and pumping losses and the speed of 1500 r.p.m. was maintained 
by partial motoring of the engine. Diagram F, shown below diagram EF, was 
taken when air only was supplied to the engine. The pressures shown accord- 
ingly can be used for estimating the approximate combustion pressures 
obtained when mixtures with hydrogen were used. 

The diagrams show that if lubricating oil passed into the combustion 
chamber at the minimum rate practicable, the rate of burning of a hydrogen- 
air mixture 5% rich or 8% weak was such that its combustion required a time 


KING AND RAND: THE HYDROGEN ENGINE. XXVII 461 


represented by nearly 15° of crank revolution and that autoignition did not 
occur even when the compression ratio was 10:1. 


V. OPERATION OF THE ENGINE AT 12:1 COMPRESSION RATIO 


The controllable operating conditions were as had been used at a C.R. of 
10, as described in Section IV, except that the temperature of the jacket 
coolant was reduced from 100° to 65°F. An unavoidable change in conditions 
occurred on raising the C.R. to 12:1, in that the top piston ring then lost its 
bearing on the cylinder wall at the three spark plug holes. Moreover a larger 
pocket into which lubricating oil could leak past the piston was formed between 
the end of the spark plug and the top piston ring than between the ends of 
the two blank plugs and the piston ring. 

The operation of the engine at a C.R. of 12 and a speed of 1800 r.p.m. in the 
conditions described was accompanied by some degree of preignition which 
occurred occasionally only and did not lead to backfiring. When however 
volumetric efficiency was increased by reducing the speed to 1500 r.p.m., 
preignition occurred occasionally with a mixture 5% rich, frequently with one 
30% rich, and developed into backfiring with a mixture 80% rich. 

The engine was then dismantled for inspection and cleaning of the combus- 
tion chamber. Patches of soot, easily wiped off, were found on the piston crown 
adjacent to the three spark plug holes. The patch adjacent to the spark plug 
hole was larger than those adjacent to the two holes occupied by blank plugs. 
Similar patches of soot were not found when the engine was dismantled 
after running at a C.R. of 10. 


Method of Avoiding Preignition by Lubricating Oil Soot and the Consequent 
Improvement in Engine Performance 

The extent of the deposits of soot on surfaces in the combustion chamber was 
found to depend on the time of running required for the accumulation of a 
quantity of lubricating oil in the pockets provided by the spark plug holes. 
Thus starting with a weak mixture and taking the observations required for 
performance data as rapidly as possible, it was found that preignition did 
not occur until a mixture richer than correct was used. When a mixture 10% 
weak was again used and the engine run for 10 to 15 min., the excess soot 
on the piston crown burned off and observations for performance data for 
the richer mixture without interruption by preignition were obtained. With 
this procedure, performance data were obtained for mixtures ranging from 70% 
weak to 80% rich and, in the absence of preignition, higher values for power 
and thermal efficiency were obtained than when experiments had been inter- 
rupted by preignition and backfiring. 

Values for I.H.P. obtained at speeds of 1200, 1500, and 1800 r.p.m. and with 
hydrogen—air mixtures ranging from 70% weak to 80% rich are given by the 
graphs of Fig. 11. 

Values for I.H.P., B.H.P., I.M.E.P., I.T.E., and opt. spark timing are 
given for a speed of 1500 with mixtures ranging from 70% weak to 80% rich, 
by the graphs of Fig. 12. 
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Fic. 11. Power developed at 12:1 C.R. on reducing the concentration of soot derived from 
pyrolysis of the vapor of excess lubricating oil and thus avoiding nuclear ignition. 
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Fic. 12. Performance data obtained in the absence of nuclear ignition when using mixtures 
70% weak to 80% rich at 1500 r.p.m. anda C.R. of 12:1. 


Performance data obtained when the mixture strength was that required 
for maximum power and that required for maximum indicated thermal 
efficiency are of especial interest and are given in Tables VI and VII respect- 
ively. Values of opt. spark timing and volumetric efficiency are also given. 

It will be seen by references to the tables that a maximum of 146 p.s.i. 
was obtained for I.M.E.P. when the engine was operated at 1500 r.p.m. and 
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TABLE VI 


MIXTURE STRENGTH FOR MAXIMUM POWER; AGED OR CLEAN EXHAUST VALVE 
HIGH SPEED CAMSHAFT, ONE COOL SPARK PLUG, COMPRESSION RATIO 12:1 





















Speed, % Vol. Spark 

r.p.m. eff. timing I.H.P. I.M.E.P. I.T.E.% 
1800 85.4 1° retard 11.9 140 37 
1500 89 1° retard 10.3 146 37 
1200 87.3 2° retard 8.0 141 36 















TABLE VII 
MIXTURE STRENGTH FOR MAXIMUM INDICATED THERMAL EFFICIENCY; 


AGED OR CLEAN EXHAUST VALVE 
HIGH SPEED CAMSHAFT, ONE COOL SPARK PLUG, COMPRESSION RATIO 12:1 











Speed, % Vol. Spark 
eff. timing 1.H.P. I.M.E.P. 1.T.E.% 











90 47 
96 47 
91 44 
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that a maximum of 47% was obtained for I.T.E. when the engine was operated 
at 1800 r.p.m. Values of 473% were not uncommon but are not shown on the 
relevant graph of Fig. 12. 

Lubricating oil always passed the piston at an abnormal rate when the 
C.R. was raised to 12. There was therefore always an appreciable concentration 
of lubricating oil vapor in the compressed hydrogen-air mixture. The pyrolysis 
of this vapor produced the finely divided carbon which remained in part as 
soot on the piston crown and required removal by the method described in 
order to avoid preignition. In these circumstances the difference in the igniting 
effect due to the exhaust valve being aged or initially clean was masked by 
the nuclear igniting effect. 

The procedure adopted to obtain performance data with a hydrogen-air 
mixture at a C.R. of 12 was made necessary by the particular design of the 
C.F.R. cylinder and piston. It could therefore be changed in a manner such 
that experiments could be carried out at compression ratios higher than 12:1. 
Thus we have, by limiting the time during which oil passed the piston, been 
able to operate for a short period at a C.R. of 16:1 with a maximum power 
mixture. 




















DISCUSSION 









Hydrogen either free or in combination with carbon is a constituent of the 
gases or liquids commonly used as engine fuels. The characteristics of the 
combustion of hydrogen only should therefore be of fundamental significance. 
It seems however to have long been accepted that hydrogen cannot be used 
as the sole fuel, because of the occurrence of premature ignition, in the operat- 
ing conditions required to obtain optimum values for power output and the 
corresponding thermal efficiency. 
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Lucke (14, 1906) states that ‘‘Gases rich in hydrogen preignite easily and 
that approximately one atmosphere should be deducted from the compression 
pressure allowable with no hydrogen for each 5% of hydrogen present”’. 
Tizard (19, p. 22, 1921) states that “Hydrogen detonates easily in an internal 
combustion engine, at quite low compression ratios, but nevertheless it has a 
high ignition temperature’. Ricardo (16, pp. 334-337, 1922), with the E35 
engine operated at 1500 r.p.m., found that it was possible to use a nearly 
correct hydrogen-air mixture at a compression ratio of 5.45:1, but “if an 
attempt was made to use a rich mixture, violent preignition occurred, accom- 
panied by firing back into the carburetor, which rendered further running 
impossible. Even with the compression ratio lowered to 3.8:1, the same 
thing occurred’’. Burstall (2, 1927) operating another E35 engine, at 1000 
r.p.m. and compression ratios of 5 and 7:1, found that ‘The experimental 
range (of mixture strength) which can be used with hydrogen is strictly limited 
by its readiness to detonate, both at high compression ratios and with rich 
mixtures, and even more by its tendency to backfire with mixtures richer than 
20% weak at a compression ratio of 5:1’’. Premature ignition, described as 
detonation, was obtained by Egerton, Smith, and Ubbelohde when they were 
using an ethyl knock testing engine (5, 1935). The effect was attributed at the 
time to nitrogen peroxide formed from the nitrogen of the air, but some years 
later to preignition (4, p. 917, 1949). 

Downs, Walsh, and Wheeler (3, p. 517, 1951), as a result of experiments with 
hydrogen as the fuel for a Ricardo E6 engine, state that ‘‘In most cases, owing 
to the extreme ease with which it ignites, a mixture richer than about 25 to 30% 
weaker than chemically correct cannot be used. If the attempt is made to run 
with a richer mixture, self-ignition occurs so early in the compression stroke 
that the mixture is eventually ignited in the induction manifold”. The state- 
ment that hydrogen—air mixtures ‘‘ignite with extreme ease’’ cannot be accep- 
ted as fact. Reference was made to the finding of King, Wallace, and Maha- 
patra that, by keeping the combustion chamber as clean as possible, the engine 
could be operated over the whole mixture range without backfiring (12). 
It was not mentioned that this occurred at a compression ratio of 10:1 and a 
speed of 900 r.p.m. The reference is also misleading in that it should have been 
to fluffy carbon that accumulated in an initially clean combustion chamber as 
described on page 269 of reference (12). Moreover it was not mentioned that 
preignition and backfiring occurred when charcoal dust was added to the 
mixture supplied to the engine. Although Ricardo’s early statement regarding 
the impossibility of using rich mixtures is echoed, they describe the premature 
ignition as ‘‘self-ignition’’ instead of by the term preignition used by Ricardo. 
Self-ignition was treated as detonation and attributed to a chain reaction 
mechanism to be described in a subsequent publication. 

Anzilotti and associates (1, 1954) describe experiments with a C.F.R. engine 
operated at 900 r.p.m. with a jacket coolant temperature of 212°F. and an 
inlet mixture temperature of 100°F. The results obtained were similar to those 
of Downs, Walsh, and Wheeler (3) in that, when the mixture strength re- 
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quired for maximum power was used, the highest compression ratio usable 
without knock was approximately 8:1. It is stated that the “knock” was 
similar to that obtained with paraffinic fuels and was attributed to the chain 
reaction mechanism devised by Lewis and Von Elbe to explain the autoignition 
of hydrogen—oxygen mixtures. In our preliminary experiments also carried out 
with a similar engine and operating conditions but with a combustion chamber 
initially clean, the knock limited compression ratio was 10:1 when a nearly 
correct mixture was used (see Fig. 1, this Part) and it was shown by subse- 
quent experiments described in Section I to have been due at 900 r.p.m. to 
the ceramic core of the spark plug having attained an igniting temperature. 

Livengood and Leary (13) obtained what appeared to be autoignition of a 
correct hydrogen—air mixture when using the compression ignition machine 
developed by C. F. Taylor and associates (18). The compression ratio was 12:1 
and after nearly adiabatic compression the mixture was held against the thick 
plate glass window which formed the cylinder head until ignition occurred 
after a few milliseconds. The window was used for photographic recording of 
the process of ignition and therefore was not cooled. In view of the experiments 
described in Section III in which a mixture 33% weak was ignited by the 
surface of the ceramic core of a spark plug at a compression ratio of 10:1, it 
may safely be assumed that the ignition observed by Livengood and Leary 
occurred because a layer of molecules on the surface of the uncooled plate glass 
window attained the required temperature. Thus it was observed by them that 
“all mixtures give off a bright yellow light during autoignition. When specimen 
inflammations were examined with a spectrograph, the presence of a strong 
sodium-D line was revealed in each instance’. Efforts were made to find and 
eliminate the source of the sodium but apparently it was not realized that 
it could be provided by the plate glass. In spite of their efforts, “‘all fuels 
including hydrogen, continued to emit the characteristic yellow light’’. 

The experiments described in this Part show clearly that optimum perfor- 
mance in terms of power and economy is limited only by preignition arising 
from the temperature attained by the ceramic core of the spark plug or by the 
exhaust valve or from nuclei formed by pyrolysis of lubricating oil vapor. 
It was shown by indicator diagrams that detonation, defined as auto or self 
ignition, did not occur in the absence of nuclei of ignition. 

On reducing the magnitude of the igniting effects mentioned it became 
possible to operate the engine at speeds rising to 1800 r.p.m., with the com- 
pression ratio fixed at 12:1. Values of 146 p.s.i. for I.M.E.P. and of 37% for 
I.T.E. were obtained when the maximum power mixture strength was used. 
The values obtained with a mixture 55% weak were 90 p.s.i. for I.M.E.P. 
and 47% for I.T.E. These values exceed those recorded in the literature of the 
subject, by a wide margin. 
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APPENDIX 


Metering of Hydrogen Supply 

Accurate measurement of rate of hydrogen supply to the engine is of prime 
importance in respect of determinations of thermal efficiency. The novel 
method used for the experiments was based on the use of two calibrated square 
edged orifices in thin plates. One was 0.1286 in. and the other 0.0784 in. 
diameter. The plates were 0.013 in. thick. The orifices were mounted midway 
in brass tubes 4 in. diameter and 4 ft. long provided with baffles and wire mesh 
screens for the elimination of pressure pulsations from the induction system 
of the engine. The double system of hydrogen supply to the engine through 
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small orifices of different diameters, which could be used singly or in parallel, 
was such that hydrogen could be supplied at rates required for the experiments 
while the velocity of flow through a metering orifice was never less than 750 ft. 
per sec. This feature in conjunction with the means taken to eliminate pressure 
fluctuations made possible the use of orifice calibrations determined in steady 
flow conditions. 

The orifices, in position in the tube mountings, were calibrated in the Meter 
Testing Department of the Consumer’s Gas Company, Toronto, under the 
direction of Mr. J. Alex Morrison, Supervisor of the Appliance Laboratory. 
A special gasometer which had been calibrated by the Ontario Department of 
Weights and Measures was used. A graph obtained accordingly for the relation 
between pressure difference and rate of hydrogen flow through the 0.1286 in. 
diameter orifice as plotted on a logarithmic scale is given by Fig. 13. The values 
given for rates of hydrogen flow are for a temperature of 60°F. and a pressure 
of 30 in. Hg. 
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Fic. 13. Calibration graph for orifice 0.1286 in. diameter. 


Accuracy of Measurement of Thermal Efficiency 

The accuracy of measurement of I.H.P. by the motoring method, which is 
in general use, is not regarded as better than +3 of 1%. It is considered that 
the rate of hydrogen consumption was made with an accuracy of +} of 1%. 
The values for I.T.E. given in this part are comparable with those obtained by 
Ricardo (16) and Burstall (2) when they used the motoring method of meas- 
uring I.H.P. Values for power and thermal efficiency are not given in more 
recent publications by others (1, 3). 


The Spark Plugs 


The spark plugs used for the experiments are illustrated by Fig. 14. Champ- 
ion C.F.R.-8 is near the top of the commercial heat scale and was used as a 
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hot variety with a spark gap of 0.025 in. Spark plugs K.L.G. 100 and R2-S, 
rated at the bottom of a commercial heat scale, were used as cool varieties, 
The spark gaps of these plugs are not readily adjustable. They were used as 
supplied with gaps of 0.018 in. approximately. It would be expected according 


















R-2-S K.L.G. 100 CFR-8 
COOL PLUG COOL PLUG HOT PLUG 


Fic. 14. Spark plugs used for the experiments. 


to the designs that the K.L.G. plug would run hotter than the R2-S but it 
was shown by experiments, Section ITI, that neither attained the temperature 
required for the preignition of hydrogen—air mixtures in the conditions of the 
experiments of this Part. 


Valve Timing 
The low and high speed camshafts were designed to provide valve timings as 


tabulated below, according to the C.F.R. handbook. The given values were 
found to be correct within 1°. 


Low speed camshaft High speed camshaft 
Inlet valve opens 10° after t.d.c. 15° before t.d.c. 
Inlet valve closes 34° after b.d.c. 50° after b.d.c. 
Exhaust valve opens 40° before b.d.c. 50° before b.d.c. 
Exhaust valve closes 15° after t.d.c. 15° after t.d.c. 
Valve lift 0.238 in. 0.312 in. 
Valve overlap 5° 30° 


When the high speed camshaft was used, the exhaust valve with a lift of 
0.312 in. was off its cooling seat for 155° of crank revolution. When the low 
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speed camshaft was used, the exhaust valve with a lift of 0.238 in. was off its 
cooling seat for 145° of crank revolution. 


Exhaust Valves 
The types of sodium cooled and uncooled silchrome exhaust valves used for 
the experiments described in this Part are illustrated by Fig. 15. 


SILCROME SODIUM COOLED 
VALVE VALVE 


Fic. 15. Types of exhaust valves used for the experiments. 
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FLUID TURBULENCE AND SUSPENDED SEDIMENT! 


By A. W. Marris? 


ABSTRACT 


A formula is derived for the concentration of fine sediment, ‘‘wash load”’, for a 


stream of infinite width flowing turbulently over a smooth bed. The theory is 
analogous to that currently employed in the problem of turbulent heat transfer. 
Sensibly exact formulae and then simple close approximations are derived for the 
total sediment load borne by the stream. 


PRINCIPAL SYMBOLS 


= concentration at distance y above bed of stream. 


ll 


Il 


I 


concentration at non-dimensional depth z. 
concentration at bed. 

mean concentration. 

depth of stream. 

reference distance above stream-bed. 

flow velocity at position y. 

mean flow velocity for stream. 

shear at position y. 

shear at bed. 

V (t0/p) = shear velocity. 

mass density of water. 

settling velocity of sediment in still water. 

Umax/V (to/p) = constant in velocity distribution law. 
von Karman’s constant in turbulent velocity distribution law. 
Nw/u, = exponent in concentration law. 

eddy diffusivity for sediment. 

eddy diffusivity for momentum. 

€,/€m- 

viscosity of water. 

u/p = kinematic viscosity of water. 

pressure at point y. 

distance measured in direction of mean flow. 
1—(y/h) = non-dimensional depth below surface. 
value of z at which concentration is at its mean value. 
u/u, = non-dimensional velocity. 

u,h/v = Reynolds number for shear velocity. 


d = boundary of laminar layer and turbulent core assuming no intermediate 
‘buffer layer’. 
6 = 1—-—(d/h). 


1Manuscript received A pril 7, 1956. 

Contribution from the Department of Civil Engineering, University of British Columbia, 
Vancouver, British Columbia. 

2Instructor in the Department of Civil Engineering. 
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= eddy diffusivity for vorticity. 
sediment velocity in direction of mean flow at P. 
sediment velocity in direction of y. 
constants of slip of water over sediment. 
a; = molecular diffusivity of sediment. 
Sc = v/a, = Schmidt number for molecular diffusion of sediment. 
S = kSec. 
a = non-dimensional boundary of laminar boundary layer. 
8 = non-dimensional boundary of “‘semiturbulent’’ buffer layer. 
L = width of stream. 
We = total sediment load carried by stream. 
W.' = W/L = sediment load per unit width of stream. 
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INTRODUCTION 









For the theoretical prediction of the concentration depth distribution of 
suspended sediment in a turbulently flowing body of water it has been cus- 
tomary to employ the formula (1, 22) 








(1) 






which expresses the ratio of the concentrations and heights y and d above the 
bed of a river of depth / in terms of these heights. The exponent @ is equal to 
Nw/u, where u, = ~/(ro/p), to being the shear at the stream-bed, p being the 
mass density of the water, V von Karman’s constant occurring in the turbulent 
velocity law (11, 13) 







u—U y 
2 ——_ = N log,* 
[ ] V/ (To ‘p) Se h ’ 
and w the settling velocity of the sediment in still water. 
The formula [1] is derived on the assumption that under steady conditions 
the upward transport of the sediment by the turbulence balances the settling 


due to gravity. One may write 


[3] wC = —e,dC/dy 










where C is the concentration of that portion of the sediment with a settling 
velocity w, and e, is the eddy diffusivity for the sediment due to the turbulence. 
e; is directly analogous to the diffusion coefficient of Fick’s law in which the 
quantity mass is transferred by microscopic molecular fluctuations, or again 
to the thermal diffusivity for gaseous heat transfer. 

Integration of equation [2] yields the relation 


*"” dy 


{ x, = —y 
| log. Ca * d €s 


in which C, is the concentration at the datum level y = d. 
The velocity of flow being a function of y, e, must also be. 
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On the statistical theory of Reynolds, Prandtl, von Karman (4, 5, 6, 12, 13, 
15) the equation for horizontal turbulent flow in the x direction may be written 


; idP_d( d fe 22 ; 
ol NE = E(u dtr Ft = BT tort Se 


pdz dy ” dy" ~ dy 
where dP/dx is the mean pressure gradient and du/dy is the velocity gradient 
perpendicular to the mean flow direction. yv is the kinematic viscosity u/p, 
u being the viscosity, and ey, is the eddy diffusivity for momentum. €),, which 
is a function of depth, is effectively an additional kinematic viscosity caused 
by the turbulence. Equation [5] is derived on the.assumption that momentum 
is transferred by the velocity fluctuations to increase the pressure gradient. 
Implicit in equation [5] is the assumption that the turbulent pressure fluctu- 
ations play no part in the over-all transfer of momentum (18, 19). 
Expressing the equations for the steady motion as 


h-y 
[6] T = T—— = 702, 


h 
where 7 is the shear stress at the point y, z = 1—(y/h) is the depth below the 
free surface in non-dimensional form, and 


[7] a Oe S: 


and defining a non-dimensional velocity 
[8] U = u/u, where u, = ~V/(r0/p) 
and a Reynolds number for friction velocity 


[9] QO = hu,/v 


one has, from equation [5], 


ro _ a aoe | | 
hp are fT (ety) dy — h’ dz (eu +v) dz ’ 


i.e. 
tn Baa PE) 
on v = 2] (41 dz 
or 
eu _ _ Qz+dU/dz 
[10] v i dU /dz 
If the velocity distribution law is given by 
[11] U= Umaxt+ N log.(1—2) 


equation [10] gives 
[12] €y/v = (Q/N)z(1—z)—1. 


It is seen that the graph of e,/v is a parabola and that e,/v is negative 
throughout intervals of length O(1/Q) at z = 0 and z = 1. The latter interval, 
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for the case of flow over a smooth bed, would be of the order of the width of 
the laminar flow boundary layer for which, of course, equation [4] does not 
hold. 

Writing ¢4/v = (Q/N)z(1—z) and proposing equality between ey and e¢; 
one obtains from [3] 


C Nwh (* dz Nef 
3 2(1—2) 


ne eS EF 


Co Qv Js2(l—2z) 4, 
where 6 = 1— (d/h), which on integration yields 


(1) 


Experiments (1, 22) have shown that [1] yields results of sensibly the correct 
form provided the exponent is 20%-50% less than the value Nw/u,, so that 
the concentration distribution is more uniform than anticipated theoretically. 
This is particularly the case for small values of @ corresponding to the very 
fine sediment. It would be expected that such would obey the law better, as 
the finer the sediment the more capable it should be of following the turbu- 
lence. It is also found that there is a tendency for the sediment content near the 
surface to exceed that indicated by the distribution law, a phenomenon at- 
tributed to random turbulence caused by the sides of the stream. Such random 
turbulence would increase e;. 

It has been stated that a statistical theory of turbulence in which momentum 
is considered conserved and regarded as the property transferred by the eddies 
requires that the fluctuating pressures that certainly exist in the turbulent 
field do not influence the mean transport of momentum. This was first pointed 
out by Taylor (18, 19), who simultaneously proved that for this reason the 
momentum transfer theory cannot hold if, for example, the motion is wholly 
two dimensional. This was verified experimentally for the case of turbulent 
flow in the annulus between rotating coaxial cylinders. He further showed 
that the only configuration under which the conditions for the momentum 
transfer theory are met is for momentum in a given direction x, say, to be 
transferred in the plane perpendicular to x by turbulent motion in which the 
flow lines parallel to x are not interfered with but remain in that direction 
throughout the motion. Such seems to be the case for flow through a straight 
tube of circular section. The theory of heat transfer, for example, based on an 
analogy with a momentum transfer theory of turbulence has been reasonably 
well ratified by experiment for this particular case. 

To overcome the difficulty for purely two-dimensional motion, Taylor (18, 
19) suggested that if the whole motion is confined to a plane yz say, then as 
in the case of two-dimensional motion of an ideal fluid, vorticity in the x 
direction is conserved, and could be taken as the property transferred by the 
turbulence to cause the pressure gradient. The theory was subsequently 
generalized to three dimensions, and a modified vorticity transfer theory set 
up, one of the main assumptions of which was that each of the components of 
vorticity was conserved and transferable in the sense of say mass in a kinetic 
theory of molecular diffusion (20, 21). 
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This modified vorticity transfer theory yields a turbulent velocity distribu- 
tion for pipe flow in as good agreement with experiment as the logarithmic 
decrement law derived from the momentum theory. The theory gives good 
results for wakes (where the turbulence may be approximately two dimensional 
in the form of a von Karman vortex street) but it breaks down in the neigh- 
borhood of solid boundaries. 

The vorticity vector being the curl of the velocity vector, the eddy diffusivity 
for vorticity involves the second differential of the velocity. It may be shown 
that on this theory the equation analogous to [5] is 


1dP d’u d’u 


du 
[13] p dz oe wt’ dy" - (ev +r) 2 


which yields, as before, by equations [6], [7], [8], [9] 


[14] ae 


vy @U/dz 
(Q/N)(1—2)*-1 


[i5] 


for the velocity distribution law [11]. 

It is seen that ey/v is negative throughout an interval adjacent to the bound- 
ary defined by 1—(N/Q)? < z < 1, an interval of width O(1/Q) and so of 
higher order than either the corresponding interval on the momentum trans- 
fer theory or the laminar boundary layer, the latter for a smooth boundary 
being of thickness O(1/(Q). 

It is seen that ey = (Q/N)—1 at the free surface z = 0 instead of zero as is 
€mu- 

The assumption of equality between e, and ey and the insertion of [15] in [3] 


now yields 
i ee Ja 
OC” ase (1—z)° 


[16] & = A; exp a where A; = exp( -;2,) ; 
Here it is assumed that 5 > O(1/0Q?). Writing (Q/N) (1—2)?—1 = (Q/N) (1—2)? 
is an approximation equivalent to neglecting v in comparison with ey. 

Whereas equation [16] certainly yields the experimentally observed finite 
concentration at the surface z = 0 (taking 6 = 0.5, 6 = 0.5, Co/Cy = exp(—}4)), 
the formula gives an infinite concentration near the wall and wholly too steep 
a concentration gradient. 

It appears then that the observed sediment distribution cannot reasonably 
be predicted on the basis of the vorticity transfer theory and one returns to 
the momentum transfer analogy. 

The author attempts to produce a more accurate result than equation [1] 
by considering, instead of a general vertical equilibrium, the transport of sedi- 
ment through an infinitesimal fluid element and integrating the ‘‘equation 


so that 
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of continuity” thus obtained. The work is roughly analogous to that of Boelter 
(2), Martinelli (10), Lion (7), Reichardt (14), and the author (8, 9) for the 
problem of turbulent heat transfer. The theory is applicable to the case of 
rectilinear flow over a smooth bed; it assumes that the river is sufficiently wide 
for the effect of the banks to be ignored. It also assumes that the sediment is 
sufficiently minute to exhibit no gravitational acceleration and to form a 
denser continuum within the continuum of the water. As in the case of radial 
heat transfer from a pipe wall to a fluid flowing turbulently through the pipe 
it is necessary to assume a linear longitudinal concentration gradient at the 
bed whereas in fact the concentration will be constant along the bed. In this 
work an uncorrected result is obtained under the assumption of the linear 
gradient and this is then used to obtain the appropriate correction. The 
method employed is that described by Seban and Shimazaki (16) for the heat 
transfer problem. 


DERIVATION OF UNCORRECTED FORMULA ON BASIS OF ANALOGY 
BETWEEN TURBULENT MOMENTUM TRANSFER AND SEDIMENT TRANSFER 


Consider a body of water flowing in steady turbulent motion in the horizon- 
tal direction x as shown in Fig. 1. The equation of motion of the water is given 
by equation [5] 


, 1dP_d du 
[5] a a (eut+v) iu |, 


Here u is the horizontal velocity at height y above the bed. 

As has been seen, equation [5] serves as a definition for ey in terms of the 
transverse velocity gradient 
[10] ew _ __Qz+dU/dz 

y dU /dz 

where U = u/V/(ro/p) = u/u, and z = 1—(y/h). 

Consider now the sediment passing through a rectangular prism of section 
ABCD as shown in Fig. 1. If u, and v, are the sediment velocities the equation 


Free Surface y=h 
v,C + que) 





D | Cc 
‘ uC coool u.C+ af uC) 
v,C 
x Bed y=0 





Fic. 1. 
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expressing no accumulation of sediment within the element is 


[17] d(u,C) 4 2U.C) = 0 
dx dy 

where 

[18] v,;C = (e,+a,) dC/dy. 


Here e, is the eddy diffusivity for the sediment, and a, is the molecular diffu- 


sivity, independent of depth. 
By [17] and [18] 


d(u;C) d dC 
[19] .* -4 (tay) se | 
and the similarity between equations [5] and [19] is at once apparent. 

Before proceeding further a relation must be assumed between the longi- 
tudinal velocity u, of the sediment and the velocity of the stream u. It will be 
assumed that there is a slip velocity between the sediment and the stream 
proportional to the stream velocity. 

Writing w—u, = j’u, one has on this assumption u, = ju where j = 1—/’. 
The motion being rectilinear and steady, du/dx = 0 and equation [19] becomes 


~~ OG. d dC 
20 = — = 
(20) I” ax a) (e+) dy | , 


In terms of the non-dimensional velocity and distance U = u/u, = u/V/ (ro0/p), 
z = 1—(y/h), equation [19] may be written 


; ae d achat | 
9) ) _ = 
2 jonu’ | 


lx dg vy dz 


where Q = hu,/v. 

It is necessary to assume now that dC/dx is independent of 2; i.e. that the 
longitudinal concentration gradient is the same for all strata parallel to the 
bed. This will not be the case in practice when the water flows over a bed of 
constant concentration. In this case dC/dx is zero for x = 1. However it is 
believed that the error introduced will only influence a small region adjacent 
to the bed. The magnitude of the correction for this will be investigated in the 
next section. 

Under this assumption equation [21] can be written 


dC poe. ee f ; 
= —jQh U dz, 
dz J€ dx Estas 0 
as dC/dz = 0 when z = 0. 
The quantity dC/dz is now eliminated from equation [22] in favor of the 
quantity (v,C),, the value of (v,C) at the bed, by means of equation [17]. 
Since adjacent to the bed there exists a boundary region of laminar flow, 


[22] 


(v,C), may be taken as —wC, where w is the fall velocity of the sediment 


in still water. 
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By equation [17] 





.d(uC) 1d = 
se ta?" 
1 1 
io)... if dC 
or (v,C), = i J 7 dz = hj ; Ue dz 
d 1 
[23] = hj x J, udz, 
liffu- 
the surface value of (v,C) being zero. 
Again if the mean flow velocity u» is defined by 
1 1 
= J suds / J as, 
0 0 
one has (v,C), = —wC,; = hjum dC/dx 
ongi- or 
ll be [24] dC ‘dx = —wC, hjtm 
‘eam . . . . . 
and substituting in equation [22] one obtains 
“o dC 2C.0/ y ) , 
25 et | I 
ymes [25] ds Ma \eo¥Gs : Udz, 
from which it is seen that slip in the assumed form has no effect on the con- 
centration gradient. Assuming a relationship e, = key between the eddy 
/p) diffusivities it is now possible to express dC/dz in terms of the fluid velocity and 
velocity gradient. 
One has 
Me Te. ; i (Hex cy 
+a, (S45) “< v TS 
- (site dz) | 1 y 
the 7 dU /dz Se 
the 1 
SQz 
1 of = se( 1-5) -35 ‘ ) 
aa ‘dz 
it 1s ; ” 
ent where Sc = v/a,, the Schmidt number for the molecular diffusion of the sedi- 
the ment, and S = kSc. 
For the case to be considered, S > 1/Q and we may write without introducing 
error (8) 
ft say" 1 (dU 
26 | i ( ). 
wi e,+a, “\dU /dz kOQz\ dz 
Equation [22] now becomes 
die ” dC_ _wl i 3 dU ) f 
a dz ~~ Rtlm\ 2 dz Uds, 
ow, so that 
ent 7 1 1a 2%... 
[28] C.-C, = Et EE | des dts, 
Rtim Jz 22 dz, Jo 
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the suffixes being inserted to indicate the order of integration. Equation [25] 
can be integrated once the appropriate velocity function is inserted in it. 
For a fully developed turbulent region the velocity law 


‘ Mm 2 
[2] “y +N log, h 


U-, 
holds, i.e. 
[11] U= Umaxt+N log.(1—2z). 


The constant N has the value 2.5 for sediment-free water, but because of the 
effect of sediment in damping the turbulence, NV increases with sediment con- 
centration up to about 3. Umax depends on the Reynolds number for friction 
velocity according to 


Umax = M+N log, Q 


where for a smooth boundary M = 5.5. 
Thus the law may be written 


[lla] U = M+N log, Q(1—2). 


From equation [lla] it is seen that 


1—(8/Q) 1—(8/Q) 
Umax— Umean = -vf log.(1—z) as / J dz, 
0 0 


assuming that the region of fully developed turbulence extends from z = 0 to 
z = 1—(6/Q), i.e. to within O(1/Q) of the boundary. One thus has 


[29] Umeon = (M—N)+N log, Q, 
[30] le = Unies = N. 


The corresponding formula for flow through a tube of circular section is 
Umax— Umean = 3N/2. 

Following von Karman the non-dimensional velocity will be taken as 
U = M+WN log.Q(1—z) for the fully developed turbulent region defined by 
the interval 0 < zs < 1—(8/Q); as U = K+J log, Q(1—z) for the buffer 
region defined by 1—(8/Q) <z<1-—(a@/Q); and as U = Q(1—z), the 
“straight’’ portion of the Poiseauille parabolic law, for the laminar adjacent-to- 
the-wall region defined by 1—(a/Q) < z < 1. See Fig. 2. Between the several 
constants a, 8, K, J, M hold the boundary relations 


[31] J=a, K = a(1—log, a), M+N log, 8 = a[1—log,(a/8)|=¥ 


expressing continuity of U and dU/dz at z = 1—(a/Q), the boundary between 
laminar and buffer layers, and continuity of U only at z = 1—(6/Q). 
The relations [31] used in connection with the experimental results of 
Nikuradse (11) and Stanton and Pannell (17) yield the numerical values 
a =J= 5, B = 30, K = 3.05, M = 5.5, N = 2.5, y = 14, for the case of a 
smooth boundary. 
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Free Surface z=0 
U=M+ Nlog,Q(I-2) 


z= |- é 
‘sciatica a ieee aa 
U=K+Jlog,Q(I-2) » 


U=Q(I-z) Bed zs| 


Fic. 2. 


Since the molecular diffusivity a, has been neglected, there is no means by 
which the sediment can be transported across the laminar boundary layer. It 
will be considered that the sediment at the bed has extended at least to the 
lower boundary of the buffer layer and as far as the integration is concerned the 
value z = 1—(a/Q) will be taken as the upper limit of z. 

It follows from equation [28] that for 0 < z < 1—(@/Q), i.e. for points in 
the region of fully developed turbulence, 


‘ we; 1—(8/Q) 
[32] C,:-C, = ha. N neath, “TM+N log, Q(i- 21)] dz,d22 


1—(a/Q) 1 1—(8/Q) 
J J tJ J [M+N log, Q(1—2:)] dz: 
1 0 


~(8/9) 22(1—22) 


+ is [K+J log, Q(1—2:)] as, | as,\ . 


—(B/Q) 


peepee these integrals in turn, 
1—(8/Q) 
J = eas, [M+N log. Q(1—21)] dzidze 


1—(8/Q) 1 
- J ———.[(M—N+N log, Q)z2—N(1—22)log.(1—22)] dze 
2 22(1—22) 


1—(8/Q) i 
= f | ar—w- log. Q) 1... 


22 


‘4 | we—-tv log.Q) (10g.(1—2)—log, s) 


+n( ys B= G/or 


2 


[33] = | ae—wew log, 0)( tog. (1 ~s) hoe, © 0 +n 
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It follows similarly that 


«1l—(a/Q) 1 + 1—(8/Q) 
[34] | oa ' { [M+N log, Q(1—2;)] dzidz2 
Ji (8/Q) 22(1—22) « 0 


#2 






= (M—N-+N log, Q) log, ®+0 ( 1), 





and that 


1—(a/Q) Ze 
[35] \ ; 7 [K+J log. Q(1—2:)] dzidze = o(3). 


J 1-(8/Q) 22(1—22) J1-a/@) 





Substituting the results [33], [34], [35] in equation [32], we have 


















ne c«¢, a af M—N-+N log.Q) (wv hadh ~<a tae, Se lel s) 
Rum 8 a 
7? r ¥ 2 a 
= (x = 4) |+0(4). 
By [31] —N log, B+a log, 5 = M—J = M-a 


so that 


[37] Gi-C, = we (M—N-+N log.Q)(N log.(1—z) +N log.-Q+M—a) 


+N (= £5)}, 


[38] 1-4 = - | 1*(log.0)*+4 2M —a—W) ) log. 0+ (M—a)(M—N)+N’ 
1 m 
2” | 


Aa’ 


Assuming that Cy = 0, i.e. zero concentration at the surface of the water, 


2 


—N((M—N)+N log,.Q) log, ei ;, 


S,, | 


equation [38] gives 


[39] 5 “= [N*(log.Q)?+N (2M —a—N) log.0+(M—a)(M—N)+4N? 2} 
so that 


o” 


N{(M—N)-+N log,Q] log, - = N° oe e 


C. rn 
C, N° (log.Q)°+N(2M—a—N) log, 0+ (M— a)(M—N)+1N' x’ 


Again since 


Cc (" a/c / f (a/Q) f (a/c ( 1 ) 
— = — a3 lz = 2 dz+O 
C1 eo Ci 0 ; 0 Cy ne Q 


it follows from [40] that 


(41) oS: N’log, QENG M—2N)+¢5N?x° 
C, N’*(log.Q)’+N(2M—a—N) log.0+(M—a)(M—N)+5N’ 9° 


[40] 
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in which are used the integrals 


J tog.(1—s) dz = —(1—z) log.(l—z)+C, 


J X5a- >» 2 —(1—2)log.(1—2) 2+. 


1 


Inserting the numerical values a = 5, N = 2.5, M = 5.5 appropriate to 
turbulent flow along a smooth bed one obtains 


1 @ 
C. [17.3+33.1 logioQ] logo Taz +625 z, 


Ci. 33.1 (logiQ)°+20.0logy~O+118  ’ 
c 14.4 loginQ+11.5 


[41a] Ci ~ 33.1 logi0Q)?+20.0(logioQ) + 11.8 ° 


It has been shown by Vanoni (22) and others that sediment in the water has 
the effect of damping out the turbulence, that is of increasing NV. 
With N = 3 instead of 2.5 equations [40] and [41] yield 





[40a] 


[17.3-447.7 logioQ] logio +9 > 5 
tia aianedagi seat eiadieea tana 
47.7 (logioQ) +20.7 (logioQ) +16.3  ’ 
20.7 logioQ+13.3 


[415] 47.7 (logioQ)?-+20.7 (logi0Q) + 16.3 * 


[400] 


Writing 

. ~$h-- feud=a,, 
2] es) = 5 
it may be shown (3) that 


a1) =, 
ng) B-[oe(oseg)] 
a(2 sin B) 10 log,\ 2 sin 10 
2 r 2 
al (2 sin x) lot “ss | toz,(2 sin =. 


2 
T 2 
$(3) = 2 ~ 2llog.2] ; 


and Tables I and II are computed for these values from equations [40] and 
[41]. Table IV gives values of log[z/(1—z)] and log(C./C,) appropriate for 


the formula 
-s) 
a 1- 


for the case 6 = 4. It is seen that the distribution ratio C,/C; and therefore @ is 
pretty well independent of Q and that the quantities are not appreciably 
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affected by variations in N caused by the damping of the turbulence by the 
sediment load. 

The graph in Fig. 3 shows, however, that log(C,/C,) is not proportional to 
log|z/(1—z)] as is required by the distribution law [1] but that @ decreases 
with z as shown in Fig. 4 for the case Q = 105. This is in agreement with the 
findings of Vanoni that the experimental value of the exponent for fine sedi- 
ment was from 20% to 40% less than predicted by the simple power law [1]. 
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Fics. 3 and 4. 


[t is also in general agreement with the results of Anderson (1) who found that 
the exponent would have to be decreased 50°%-60% in order to fit the results. 
Sediment distribution in practice being so subject to extraneous disturbances, 
the results are naturally scattered so that it is extremely difficult to extract 
from them any conclusive information concerning the depth variation of the 
exponent #. What can be said, however, is that the distributions are generally 
more uniform than predicted by [1], and that this could well correspond to an 
“‘annulment’’ caused by @ decreasing with depth below surface (increasing 2). 
The highest values of @ occur just below the surface and Vanoni (22) repeatedly 
finds that the concentration is unduly high in this region. 
From equations [41a], [41)] it is seen that as Q becomes large 


Ca 1 l 
{41c] oo = 9 anne ’ 
Ci log.Q 2.3026 logioQ 
This approximation gives results about 3% low throughout the range Q = 10° 


to Q = 2X 10°. See Table III. 
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TABLE I 
CasE N = 2.5 


y the 


al to 


>ases C2/Ci 


" 
to 
x 
= 
So 

= 


1 the Q=10? Q=5xX10? Q=10! Q=2X10! Q=5X104 Q=7.5X10! Q=105 Q 
sedi- ; 7 ‘ _ 
r (1). 0. 


0. 
0. 
0. 
0. 
0. 
0. 
1.0 





.047 . 04 .043(5) 
.067 R 062(5) 
092 ; .086(5) 
215 ; . 202 
278 “ . 262 
424 0.406 .400 
.631 ; .595 

0 ‘ 0 


Cm/Ci 


coco 
® to 
a2 es 
uO oe 


a 
be 


kroooooocco 


eoooceo 
ow 





Q=2xX104 Q=5xX104 Q=7.5X104 Q=105 Q=2X105 


Q=103 Q=5X103 Q=104 


0.148 0.120 0.110 0.103 0.095 0.091 0.089 0.084 





TABLE II 
CasE N =3 





2 Ci 


Q = 105 Q = 2x105 


Q = 10° 


0 0 0 
.055 0.043(5) 0.041 
.079 0 .062(5) 0.059 
. 109 0.086(5) 0.082 
. 204 0.202 0.191 
.327 0.262 0.249 
497 0.400 0.378 
. 740 0.595 0.562 
0 1.0 1.0 

Cn/ Cs 


Q = 108 Q = 108 Q = 105 = 2x10 


0.149 0.110 0.089 0.084 


TABLE III 


COMPARISON OF EXACT VALUES OF C,,/C; FOR N=2.5 WITH VALUES OF 1/log.Q 


Q 108 5X 108 10* 5X10* 108 2X 10° 


0.148 0.120 0.110 0.095 0.089 0.084 


0.145 0.118 0.108 0.092 0.087 0.082 
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TABLE IV 








log Cs 
C 


Q = 103 Q = 104 Q = 105 Q = 2x105 





z log - 
8 5 


~ 
o 
g 


o 


























Case N = 2.5 
0.382 —0.209 —0.158 —0.157 —0.157 —0.158 
0.500 0 0 0 0 0 
0.618 +0.209 0.140 0.142 0.141 0.140 
0.900 +0.954 0.504 0.506 0.510 0.510 
0.950 1.279 0.615 0.617 0.622 0.625 
0.990 1.996 0.796 0.799 0.806 0.805 
0.999 3.000 0.968 0.972 0.979 0.979 

Case N = 3 
0.382 —0.209 —0.155 —0.157 —0.157 —0.158 
0.500 0 0 0 0 0 
0.618 +0.209 0.140 0.142 0.141 0.140 
0.900 +0.954 0.501 0.507 0.510 0.510 
0.950 1.279 0.613 0.617 0.622 0.625 
0.990 1.996 0.794 0.799 0.806 0.807 
0.999 3.000 — 0.971 0.979 0.979 





CORRECTION FOR UNIFORM LONGITUDINAL SEDIMENT DISTRIBUTION 
AT BED 
For a linear longitudinal concentration gradient at the boundary, by 
equations [38] and [41], 
C,—C, 


1 co 
[44] —™ Ri—R: loge 7; Med 


This result will now be approximately corrected so as to apply to the case of 
constant longitudinal concentration at the boundary. The method is that sug- 
gested by Seban and Shimazaki (16) for the heat transfer problem but used 
by these authors to obtain results appropriate for a computer rather than to 
obtain an explicit equation. 

It will be assumed that the concentration (C;—C,)/(Ci—C,,) is invariable 
in the direction of flow, that is, a function of z only, so that 


d A-C, 1 (sc: aC.) C.-C, (2 Cn) 7 


2 


n 


0 


~ @eenHe, CHa ee Ht 


If the wall concentration is constant so that dC,/dx = 0, equation [45] 
gives 


dx dx 


dC, _ Gi—C. dCn 
[46] & "C=C, &’ 


so that in place of [22] one has 


dC _ a pee * Gi—C, dC r 
[47] dz kQh €stdsJo Ci—C, dx Uds, 
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v__dCn Gi-C. 
esta, dx Jo U Ci—Cn ds, 





imnele [48] = —kQh 
and in place of [23] 


1 
(v,C), = —wC, = ie | was 


ia [49] = hk oon = fu Ge 


Now (Ci—C,)/(Ci—C,) is near unity for z < 0.9 say, i.e. through a large 
proportion of the range of integration, so that [w(Ci1—C.)/(Ci—Cn)]mean May 
be taken as Um| (Ci— C.)/(Ci- C. sa) tacit Or Um, 


pain i.e. n= fu i = as / f dz = fus ———— de, 


so that in place of equation ae 
, [50] dC,,/dx = —wC,/hkum. 


) By equations [48] and [50] as before, 


. _ _wli f 1dU a—-C, 
ON [51] C1 C: - Rum z 22 dz "U Ci- C ds 


so that finally, substituting from 4 a and approximating as before, 








, by . 1—(8/Q) 
[52] CG—C, = - 4 f mee J [M+N log.Q(1—2)] 
Rum 2 x 
x es Ro log, =— 7 Rb #) dz,d22 
se of 1—(a/Q) 
J fe M+N log.Q(1- 
sug- ? 1-(8/0) Z2(1—z2) [ ~ og Q( z)] 
used 1 2" 
ot x | Ri—R: log, mes =; | dzidz. ¢° 
able Whereas the integral 
=. 3" 
log.(1—z) z. mt 
1 
can be evaluated in terms of logs and the series )~ oe > Be the result is 
[45] a, . 1” 


very o-eeaniamaaass accordingly an approximation will be made using instead 


of Ss , the value of this series at z = }. 
: SQ 2") S (1-2)" r 
Since O(s) = >> 2 z A +log.s log.(1—2) = @(1) = |, 
1 1 
2 


YL 


= jp B(log.2)°, 
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to 





so that writing oe (3) ao R, 
1 n 
one has 
1—(8/Q) 1 Zs ; ( 1 - me ) : 








a 1—(8/Q) 1 ie " ’ 1 
j 1 “(+N log. (1—21))( Ri’ Re loge =~) de 
ewe 9 — 


22(1—22) 

sf tog. —z)—log, £\_s| » =4 
co 1 5 n 

2) > a») 

i 


=s)_ pa (3) 
n 2D n° +0 Q 


| S, = (M—N+N log.Q)Ri’ —(M—2N-+N log.Q)R:, 
[54] So = NR,'+(M—2N+N log.Q)R», 
| Ss = NR. 


ll 


[53] 


+5) (log, (1—z) ]"log,z+2 log, (1— 


where 


Putting 8 = aandz = 1—(8/Q), the second integral in equation [52] becomes 
simply 


[55] a log. © +0 (4) 


so that by [50] 


a 2 « n 
[56] C.-C, = wey SiN log-0+M—a]+N Si —NS, log, a —NS:>> =3 


-_ . =. (l—z)” = (1—3z)" 
+N S| tog, (1—z)]"log.z+2 log,(1—z) >> U=% -—=—% a=s" |} 
1 ) 
in which is used the boundary relation 


—N log, B+a log, - = M-—a. 


RID 


By equation [36] 


[57] Ri = GG (M—NN log.Q) (N log.Q+M—a)+jN’r’], 
m 1 \m 

[58] R; = oo ate (M—N-+N log.Q), 

” wc, _N* 

[59] s Rum Ci— Crh . 
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so that 


2 
Ry = Ri- Es 9 (log.2) ‘le, 


wc, 1 


” in =| w- N-+N log.Q) (N log-Q+M—a) 


- w(S43 + (log.2) ») | 
* we; 1 2 r 1 2 ry y 
[60] = 1+F ee C |» (= 4409.2) ) +N(M—-2N+N 108.0) | 


by equations [39] and [41]. 
It follows that 


S, = (M—N-+N log.Q) Ry’ — (M—2N-+N log.Q) Rz 


oe we, 1 ») | 
= (M—N+N log.)| 1426 C.-C, n( E+ +5 (log.2) . 
So = NR1'+(M-—-2N+N log.Q) ie 


“ wi __1 of x 1 ) 
[62] wy 9 abel» (H+ 4 og.2) 


\ 
+(M-+N log.Q)(M—2N+N 10¢.0) | 9 
and S3 = NR>2 


3) = ae N-+N log.Q). 


Substituting the results [61], [62], and [63] in equation [56] and writing from 
[36] 


(C1—Cz)une = wel (ay — —N+N log.Q)(N log.(1—z) +N log-Q+M—a) 


one has finally 


[64] (Ci—C-) cor = (Ci— C:)une+a 2 ( a Gt) 


aay rien 
ta,{ 26 (Ci— Calero 6 


- (=. cbs 9 (loge2 2") 


= (M+N log.Q)(M—2N+N log,Q)N?. 
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Taking M = 5.5, N = 2.5, the correction factors [65] and [66] become 
[65a] a; = 9.475, 

[66a] a2 = 17.2+216 logi0Q+207 (logioQ)’, 


and equation [64] may be written 
[67] (&) = (=) .- a(;”-)( 1-C,/C, ) 
Ci cor Ci une Rum i- Cn /Cy une 


w \? 7S Ae * < o 
-a( 2) '7 (. —C,,/C1) oy “ip n 


where (Ci)cor = (Ci)unc, i.e. the same boundary concentration is considered in 
both cases, and where by [39] and [41] 


wi 38 


; we e 7 
[68] ae (.v'te, Q+N(M—2N)+N’ 


[68a] on (14.4 logyQ+11.5)~* 

1 
for the values of M and N considered. 

Tables V, VI, and VII now give the corrected values of the concentration 
ratio C,/C,. 

It is seen that the boundary condition Cy) = 0 has had to be abandoned 
because of the approximation so that the calculated value of the correction 
must be too large near the surface. For values of Q > 10‘ the correction is 
entirely insignificant. The change in concentration between the bed and the 
fluid is localized in the laminar adjacent-to-bed region to leave the concentration 
gradient throughout the turbulent core sensibly unchanged from its almost 
constant value. 


TABLE V 
CasE N = 2,5, Q = 10%, 


W 
i @ 2.71X10-3, az = 1286, 


Cn 
G 
(c) “(¢ a. ail wor a). ere ta 
cor unc 6 1 n? 


= 0.148 





(24) (2. (2, (2. 





0 1.645 0 1 —0.021 
0.382 1.219 0.073 0.927 0.057 
0.500 1.062 0.105 0.895 0.091 
0.618 0.890 0.145 0.855 0.132 
0.900 0.5 0.335 0.665 0.327 
0.950 — 0.433 0.567 0.427 
0.990 — 0.657 0.343 0.653 
0.999 — 0.975 0.025 0.975 

— 1.000 1.000 
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TABLE VI 
Case N = 2.5, Q = 10! 


Cm = 09.110, W_ = 1.59x10-%, a: = 3312, 
Ci ktm 


C e = ™ ¢ = x? @ gn 
= =(— —0.1 10>? (1—-=* —0.9410-? (— -—z— 
a). oe : - ( Cy unc _ f: 1 n* 


x 2m c € G 
mie c. 1-(& Cc: 
(j 1 n? ).. (2)... (=) 


1.645 0 —0.017 
1.219 0.055 945 .042 
921 .067 


1.062 0.07 
0.890 0. 890(5) .099 
0.5 0.2: 747 .247 
0.: 673 .326 
0.497 503 .496 
.259 741 


0.7 
1.000 .000 











coossocs 


TABLE VII 
CasE N = 2.5, Q = 105, 


Cn — 9,089, “ = 1.065x10-3, a2 = 5175, 
Ci kum 


(2) - (2) —0.11x10-*(1— —0.60(6) x10"*( = —2 22) 
cor unc C 6 1 n? 


C1 


1/ une 


CH) @. +O. © 
6 i Ci unc Ci unc C, cor 


0 1.645 1.000 
0.382 1.219 0.043(5) 0.956(5) 
0.500 1.062 0.062(5) 0.937 (5) 
0.618 0.890 0.086(5) 0.913(5) 
0.900 0.500 0.202 0.798 
0.950 — 0.262 0.738 
0. 
0. 
0. 


C 








0.990 — 0.400 600 
0.999 — 0.595 405 
1.000 0.000 1.000 000 


TOTAL LOAD CARRIED BY TURBULENT STREAM 


It remains to consider the total wash load carried by the stream, i.e. to 
obtain a maximum value for the sediment deposited if the turbulence were 
to be stopped by flow into a reservoir for example. 

Consider a sheet of the stream contained between the planes z and z+dz. 
Then if Z (assumed large in comparison with the depth) is the width of the 
stream, the weight of sediment crossing a given section of the sheet per unit 
time is u,C,Ldz. This is the weight of sediment that in one second would be 
set free to fall back to the bed if the turbulence between the planes were to be 
suddenly stopped. 
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Assuming u, = ju where wu is the velocity of the stream at z, one has for the 
total weight of sediment transported in unit time across a plane x = const. 
perpendicular to the mean flow direction 


1—(a/Q 1—(a/Q) 
[69] Wo= j u,C,Ldz = f ju,UC,Ldz 
0 0 
wa 1—(8/Q) 1 
[70] = ju, J UCdz+0O ( 0 


—(B/Q) 


(M+N log, o2))(&) dz+O ( 1) 


where C,/C; is given by equation [40]. 
This integral again requires the integral 


oo oa” 
fos. (1—2) m2 mie 
but as the term involved is small the —— is once again simplified by taking 


2)’. 


vl 
[71] = ju,LC, J 


0 


a o” co (4)”" 
>> =; at the constant value > 44 = 
3 


1” 1 


2 —3(log. 
With this approximation the integral [71] gives 
[72] Wo = jusLC J ___N*(log.Q)'+N*[2M—3N+N (iy — 3 (log.2)") Jlog.Q_ 
a1 a = IM) NF (log.0)?-+N (2M —a—N)log-0+(M—a)(M—N)+iN' x? 
‘ N(M-N)[M—2N+N (is9" — }(log.2)’)] m 
N®(log.0) +N (2M—a—N)log.0+(M—a)(M—N)+iN'#') * 


and inserting the numerical values N = 2.5, M = 5.5, a = 5, 


We _ 82.8(logi0Q)?+71.4 logiw0Q+14.7 
ju,LC, 33.1 (logioQ)*+20.0 logioQ+11.8 


and for the higher value of VN, N = 3, 


We _ 148. 2(logi0Q)” +77 .7 logioQ+9.37 
jurLC, 47, 7 (log10Q)? +20.7 logi00+16.3 


For large values of Q, Wo/ju,LC,; tends simply to the value JN, and is 
independent of Q. The alteration of N due to the damping of the turbulence 
by the sediment load now becomes an important factor. 

The formulae [72a], [726] give the load set free per second on stopping the 
turbulence in terms of the concentration C, at the bed 
WR Assuming no slip between the sediment and the water, 7 = 1, so that if 
W.’ = W@/L the mass deposited per unit width of stream 


Wa’ _ 82. 8 (logi0Q)” +71.4 logioQ+ 14. 7 





[72a] 


[725] 


| ; ditt ee ee oe Ok 
78s} Civ, 33.1 (logiQ)?+20.0 logioQ+11.8 tor NW = 3.5 
and 

r 9 ( 
[73b] We’ _ 148.2(logwQ) +77.7 logu0+0.37 for N = 3. 


Cu, 47. 7 (logioQ)°+20.7 logi00 + 16.3 
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In terms of the mean concentration for the case of no slip, equations [72] 
and [41] give 














74] Wes iu *(log.Q)’-+N*[2M—3N + (i34" — 3 (log.2)*)N]log.Q 
l ike. N’log-O+N(M—2N)+iN' wr? 
4, N(M—N)[M-2N+ (sm 3 (log.2)*)N]\ 
N*iog.O+N(M—2N)+1N' 7, 
i.e. 
We’ _ 82.8(logiQ)°+71.4 logi0Q+14.7 —* 
[74a] Cir 14.4 logwQ+11.5 for N= 3.8, 
We’ _ 143.2(logi0Q)*+77.7 logi0Q+9.37 a) 
i Cs,” 20.7 logiv0-+13.3 et ee 
It is seen that Wg’/u,;C, approaches the value N log, Q when Q becomes 
large. 


C, can be determined from readings across the flow profile, or taken as the 
concentration at a point z,, where, by equations [37] and [38], 











: _ gg a es | 

ey alicia, [ N*log.0+N(M—N) 

oi os” logwO+115 | _— 
[75a] = J exp] 14.4 logunO-+7.5 for N = 2.5, 
“ tec a logw0-+ 13:3 ] _ 
750] = 1 exp| 20.7 logiwO-+7.5 for N = 3. 


Table VIII gives values of W9’/Ciu, and W’/C,,u, and also z,, the value 
of z at which the concentration is at its mean, for values of Q from 10° to 10°. 

It is seen that throughout the range of Q a satisfactory approximation is 
obtained by setting 


[76] W'/Ciu, = N 


TABLE VIII 

















Wa’ Wa’ N] a 

Q Ciu Cnt N log. fies 
Case N = 2.5 

108 2.63 17.8 17.3 0.660 

104 2.62 23.5 23.0 0.654 

105 2.60 29.2 28.8 0.650 

108 2.58 34.9 34.6 0.647 
Case N =3 

108 3.02 20.3 20.7 0.662 

104 3.03 27.3 27.6 0.655 

105 3.03 34.0 34.5 0.651 

108 3.03 41.0 41.4 0.648 
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where N is von Karman’s constant occurring in the original turbulent velocity 
distribution law. 

More apparent still is the closeness of Wa’/C,,u, to N log, Q throughout this 
extended range of Q. 

Also apparent is the invariance of z,,, the point at which the concentration 
takes on its mean value. The value z, = 0.65 will give results accurate to 1% 
over the range of Q. 

The concentration varying slowly over the region about z = 0.65, this seems 
a satisfactory method of determining Cp. 

Once C,, is determined experimentally the formula 



































[77] W,’ = Cy uU, log. Q 


is a simple expression for determining the deposit rate per unit stream width 
on stopping the turbulence. The friction velocity u, is directly obtainable from 
the hydraulic grade line and Q is merely the Reynolds number constructed 
with friction velocity u,. 


SUMMARY 


An equation has been derived for the sediment concentration ratio as a 
function of depth on the momentum transfer analogy for the case of infinitely 
fine sediment suspended in a river flowing turbulently over a smooth bed. The 
equation is of the form 


l oO 2” 
(a+b log.Q) log, i -+C>> 2 
—3 


[78] ae d(log.Q)°+e log.O+f 
where the constants a, b, c,d, e, f are given by the constants defining the velocity 
distribution. The velocity distribution considered consists of a laminar flow 
layer adjacent to the river bed, a semiturbulent buffer layer, and the major 
body of the fluid throughout which the turbulence is fully developed and in 
which the flow velocity is given by the von Karman velocity defect law. 

Results from the above formula are compared with those from the power 
law 


1 C, _& : a 
(1 & (st x 5 


obtained on the simple assumption of equilibrium between the upward trans- 
port of sediment by turbulence and its settling due to gravity. 

Full agreement between the two equations can only be obtained by allowing 
the exponent @ of equation [1] to vary, but approximate over-all agreement is 
obtained by taking @ at about half its value Nw/u,. Both the author’s equation 
[78] and equation [1] with the lower index imply the more uniform concentra- 
tion—depth distribution that has been found experimentally by Vanoni (22 
and Anderson (1). 
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The equation yields a simple closely approximate formula for the ratio of the 
mean concentration to the concentration at the bed, 
Co 1 1 


[41c] Ci log.Q ~ 2.3026 logi0Q* 








In the author’s original derivation it was necessary to assume that near the 
boundary, the sediment concentration, instead of being constant in the flow 
direction, had to vary linearly with distance. Using this result a first correction 
was obtained towards the constant boundary concentration case. The correc- 
tion was found to be entirely negligible; the longitudinal concentration dis- 
tribution in the adjacent-to-bed region does not appear to affect the distribu- 
tion in the main body of the water. 

It was found possible to obtain an expression for the total sediment load 
carried by the flow, that is, the sediment that would be freed to fall on stopping 
the turbulence. Without knowledge of the mechanism of “‘pick up’”’ from the 
bed (and of the transport through the laminar boundary layer) it is possible 
only to express this quantity in terms of either the concentration at the bed or 
the mean concentration. However, the simple result 


(77 Wa’ = Cy u, log. Q 


for the load per unit stream width would seem to be satisfactory for practical 
application, particularly as the mean concentration is shown to occur at the 
sensibly invariable depth 0.654, where / is the total depth of the river. 
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CORRECTION 


Vol. 33, p. 353. In the first line of Item 4 under Procedure (Reagents) “0.75 gm. sali- 
cylic acid” should read “0.075 gm. salicylic acid’’. 








